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ABSTRACT 


The  performance  of  coherent  and  noncoherent  RAKE  receivers  over  a  fading 
channel  in  the  presence  of  pulse-noise  interference  and  additive  white  Gaussian  noise  is 
analyzed.  Coherent  RAKE  receivers  require  a  pilot  tone  for  coherent  demodulation.  When 
a  first  order  phase-lock-loop  is  used  to  recover  a  pilot  tone  in  the  presence  of  additive 
white  Gaussian  noise,  phase  distortions  at  the  phase-lock-loop  output  result.  These  phase 
distortions  produce  an  irreducible  phase  noise  error  floor  for  soft  decision  Viterbi  decod¬ 
ing.  Both  coherent  and  noncoherent  RAKE  receivers  optimized  for  additive  white  Gaus¬ 
sian  noise  perform  poorly  when  pulse-noise  interference  is  present.  When  soft  decision 
convolutional  coding  is  considered,  the  performance  degrades  as  the  duty  cycle  of  the 
pulse-noise  interference  signal  decreases.  The  reverse  is  true  for  hard  decision  Viterbi 
decoding,  since  fewer  bits  experience  interference  and  bit  errors  with  high  noise  variance 
cannot  dominate  the  decision  statistics.  The  soft  decision  RAKE  receiver  optimized  for 
pulse-noise  interference  and  additive  white  Gaussian  noise  performed  the  best  for  both 
the  coherent  and  noncoherent  RAKE  receivers.  This  receiver  scales  the  received  signal  by 
the  inverse  of  the  variance  on  a  bit-by-bit  basis  to  minimize  the  effect  of  pulse-noise  in¬ 
terference.  The  efficacy  of  this  technique  is  demonstrated  by  analytical  results,  which  re¬ 
veal  that  this  receiver  reduces  the  probability  of  bit  error  down  to  the  irreducible  phase 
noise  error  floor  when  pulse-noise  interference  is  present.  This  demonstrates  how  impor¬ 
tant  it  is  to  design  the  receiver  for  the  intended  operational  environment. 
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EXECUTIVE  SUMMARY 


The  performance  of  coherent  and  noncoherent  RAKE  receivers  for  a  signal  trans¬ 
mitted  over  a  fading  channel  in  the  presence  of  pulse-noise  interference  and  additive 
white  Gaussian  noise  is  analyzed.  Unlike  commercial  systems,  which  operate  under  strict 
FCC  regulations,  military  communication  systems  are  often  required  to  operate  with  high 
levels  of  noise  and  interference  in  addition  to  frequency-selective  fading.  Reliable  com¬ 
munication  over  a  frequency-selective  fading  channel  can  be  achieved  by  either  using  a 
RAKE  receiver  or  orthogonal  frequency-division  multiplexing  (OFDM).  This  dissertation 
analyzes  the  performance  of  a  RAKE  receiver  with  convolutional  coding  with  both  hard 
and  soft  decision  decoding.  Coherent  RAKE  receivers  require  a  pilot  tone  for  demodula¬ 
tion.  The  use  of  a  first  order  phase-lock-loop  to  recover  a  pilot  tone  with  additive  white 
Gaussian  noise  causes  phase  distortions  at  the  phase-lock-loop  output,  which  produce  an 
irreducible  phase  noise  error  floor  for  soft  decision  Viterbi  decoding.  These  phase  distor¬ 
tions  are  negligible  when  the  pilot  tone  signal-to-noise  ratio  is  20  dB  or  greater.  If  the  pi¬ 
lot  tone  is  received  with  a  lower  signal-to-noise  ratio,  reliable  communication  may  not  be 
possible.  Analytical  results  for  soft  decision  Viterbi  decoding  with  a  10-dB  pilot  tone  sig- 
nal-to-noise  ratio  predicts  a  catastrophic  50  percent  error  rate.  This  demonstrates  how 
critical  it  is  for  a  coherent  communication  system  to  have  a  good  phase  reference. 

Commercial  RAKE  receivers,  such  as  the  ones  used  in  IS95  cellular  phones,  are 
optimized  for  additive  white  Gaussian  noise.  This  dissertation  demonstrates  that  both  co¬ 
herent  and  noncoherent  RAKE  receivers  optimized  for  additive  white  Gaussian  noise  per¬ 
form  poorly  when  pulse-noise  interference  is  present.  When  soft  decision  Viterbi  decod¬ 
ing  is  considered,  the  performance  degrades  as  the  duty  cycle  of  the  pulse-noise  interfer¬ 
ence  signal  decreases.  The  reverse  is  true  for  hard  decision  Viterbi  decoding,  since  fewer 
bits  experience  interference  and  bit  errors  with  high  noise  variance  cannot  dominate  the 
decision  statistics.  Soft  decision  receivers  provide  a  decision  if  the  received  bit  is  a  logi¬ 
cal  “1”  or  a  logical  “0”  as  well  as  an  estimate  of  the  of  the  receiver  confidence  in  the  bit 
decision.  If  the  pulse-noise  interference  causes  a  bit  error  and  the  soft  decision  receiver 
gives  the  bit  decision  a  high  confidence  estimate,  then  soft  decision  Viterbi  decoding  will 
produce  significant  errors.  Hard  decision  Viterbi  decoding,  on  the  other  hand,  provides 
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some  immunity  to  both  pulse-noise  interferenee  and  phase  noise.  When  the  pulse-noise 
interferenee  had  a  one  pereent  duty  eycle,  the  eoherent  RAKE  reeeiver  with  hard  decision 
Viterbi  decoding  produced  an  error  rate  of  approximately  10 or  lower  when  the  pilot 
tone  signal-to-noise  ratio  was  20  dB.  If  the  pilot  tone  is  received  with  only  a  10-dB  sig- 
nal-to-noise  ratio,  then  a  coherent  RAKE  receiver  with  hard  decision  Viterbi  decoding 
could  still  provide  reliable  communication.  An  error  rate  of  10  *’  can  be  obtained  when 
the  ratio  of  the  received  signal  power-to-pulse-noise  interference  power  is  30  dB  or 
greater. 

Noncoherent  RAKE  receivers  do  not  track  the  phase  of  the  received  signal,  hence 
they  are  not  affected  by  phase  noise.  Using  soft  decision  Viterbi  decoding,  the  noncoher¬ 
ent  RAKE  receiver  can  perform  below  the  coherent  RAKE  receiver  phase  noise  error  floor 
when  the  ratio  of  the  received  signal  power-to-pulse-noise  interference  power  is  large. 
However,  the  coherent  RAKE  receiver  for  small  values  of  received  signal  power-to-pulse- 
noise  interference  power  performs  significantly  better  than  the  noncoherent  RAKE  re¬ 
ceiver  when  the  effects  of  phase  noise  are  negligible.  This  is  especially  true  when  hard 
decision  Viterbi  decoding  is  considered.  For  a  pulse-noise  duty  cycle  of  one  percent  and 
approximately  the  same  received  signal  power  and  pulse-noise  interference  power,  the 
probability  of  bit  error  for  a  noncoherent  RAKE  receiver  with  hard  decision  Viterbi  de¬ 
coding  was  approximately  10 .  This  is  four  orders  of  magnitude  worse  than  the  coherent 
RAKE  receiver  with  hard  decision  Viterbi  decoding. 

A  RAKE  receiver  optimized  for  pulse-noise  interference  and  additive  white  Gaus¬ 
sian  noise  performed  the  best  for  both  the  coherent  and  noncoherent  cases.  This  receiver 
scales  the  received  signal  by  the  inverse  of  the  variance  on  a  bit-by  bit  basis  to  minimize 
the  effect  of  pulse-noise  interference.  The  efficacy  of  this  technique  is  demonstrated  by 
analytical  results,  which  reveal  that  this  receiver  reduces  the  probability  of  bit  error  down 
to  the  irreducible  phase  noise  error  floor  when  pulse-noise  interference  is  present.  This 
demonstrates  how  important  it  is  to  design  the  receiver  for  the  intended  operational  envi¬ 
ronment. 
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I.  INTRODUCTION 


Wireless  eommunications  over  multipath  fading  channels  are  increasingly  impor¬ 
tant  for  modem  digital  communications  systems.  This  is  especially  tme  for  military 
communication  systems,  which  may  also  have  to  operate  in  high  levels  of  noise  and  inter¬ 
ference.  There  are  many  examples  in  modern  digital  communication  systems  of  wireless 
transmission  over  fading  channels,  such  as  cellular  telephones,  wireless  local  area  net¬ 
works,  and  the  physical  layer  of  a  mobile  ad  hoc  network,  such  as  a  network  between  two 
ships.  Furthermore,  with  the  demand  for  even  higher  data  bit  rates  to  support  ever  more 
sophisticated  data  applications,  the  fading  channel  is  more  likely  to  be  frequency- 
selective.  Consequently,  modern  digital  communication  systems  have  to  be  designed  to 
operate  over  frequency-selective  channels  in  order  to  avoid  significant  degradation  in 
performance. 

One  technique  that  can  be  used  to  eliminate  the  degradation  that  results  from  a 
frequency-selective  channel  is  orthogonal  frequency-division  multiplexing  (OFDM) 
where  the  data  bits  are  transmitted  on  a  number  of  orthogonal  carrier  frequencies.  OFDM 
is  a  special  case  of  multicarrier,  as  opposed  to  single  carrier,  modulation  where  the  use  of 
orthogonal  carrier  frequencies  minimizes  overall  signal  bandwidth.  If  there  are  N  or¬ 
thogonal  carrier  frequencies,  then  one  out  of  every  N  data  symbols  is  transmitted  on  the 
same  carrier  frequency,  and  the  symbol  rate  for  an  individual  carrier  is  reduced  by  a  fac¬ 
tor  of  N.  As  a  result,  the  bandwidth  required  for  each  carrier  is  reduced  by  a  factor  of  N  as 
compared  with  the  bandwidth  of  the  signal  when  only  a  single  carrier  frequency  is  used. 
By  choosing  A  large  enough,  the  channel  for  each  carrier  can  be  made  frequency- 
nonselective,  or  fiat,  even  though  the  channel  is  frequency-selective  when  only  a  single 
carrier  is  used.  OFDM  essentially  works  by  converting  a  single  signal  into  multiple  sig¬ 
nals,  where  each  of  the  multiple  signals  is  frequency-nonselective  over  the  same  channel 
for  which  the  original  single  carrier  signal  is  frequency-selective. 

Another  technique  for  minimizing  the  performance  degradation  caused  by  fre¬ 
quency-selective  channels  involves  the  modification  of  the  receiver  rather  than  the  trans¬ 
mitter.  One  type  of  receiver  designed  to  demodulate  signals  that  have  been  transmihed 
over  a  frequency-selective  channel  is  the  RAKE  receiver.  In  this  dissertation  the  perform- 
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ance  of  a  RAKE  receiver  when  the  signal  is  transmitted  over  a  frequency-selective  chan¬ 
nel  is  analyzed.  Since  military  systems  face  even  more  complicated  issues  such  as  inter¬ 
ference,  the  effect  of  pulse-noise  interference  as  well  as  additive  white  Gaussian  noise 
(AWGN)  is  also  analyzed.  For  coherently  detected  signals,  the  effect  of  phase  noise  on 
the  system  performance  is  analyzed.  A  block  diagram  of  the  assumed  transmitter  is 
shown  in  Figures  1.1. 


Figure  1 . 1  Transmitter  Diagram. 

As  can  be  seen  in  Figure  1.1,  the  data  is  assumed  to  be  convolutionally  encoded,  inter¬ 
leaved  and  modulated  as  either  a  Binary  Phase  Shift  Keying  (BPSK)  waveform  or  fre¬ 
quency  shift  keying  (FSK)  waveform.  The  signal  is  assumed  to  be  transmitted  over  a 
slow,  frequency-selective  fading  channel  with  both  pulse-noise  interference  and  AWGN. 
Without  loss  of  generality,  in  order  to  evaluate  the  effect  of  channel  fading  and  pulse- 
noise  interference  on  the  data  signal  alone  as  well  as  on  both  the  data  signal  and  the  pilot 
tone,  the  pilot  tone  is  assumed  to  be  transmitted  on  a  separate,  independent  frequency 
from  the  data  waveform  carrier  frequency. 

A  block  diagram  of  the  receiver  is  shown  in  Figure  1.2. 


Figure  1.2  Receiver  Diagram. 


A  RAKE  receiver  is  used  to  minimize  the  effect  of  the  frequency-selective  fading  chan¬ 
nel,  and  the  maximum-likelihood  RAKE  receiver  for  AWGN  and  pulse-noise  interference 
is  analyzed  and  compared  with  a  RAKE  receiver  designed  for  AWGN  only.  De¬ 
interleaving  the  signal  randomizes  the  errors  prior  to  decoding.  Both  hard  decision  and 
soft  decision  Viterbi  decoding  are  analyzed.  In  addition,  phase  noise  effects  are  also  con¬ 
sidered  for  coherent  systems. 
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The  effects  of  phase  estimation  errors  in  coherent  RAKE  receivers  have  been  ex¬ 
amined  in  the  literature.  In  [1]  and  [2]  a  RAKE  receiver  operating  at  a  low  signal-to-noise 
ratio  (SNR)  was  analyzed.  At  low  SNR  the  RAKE  receiver  has  problems  estimating  the 
channel  tap  weights,  and  at  any  given  time,  each  finger  will  have  a  different  random 
phase  error.  A  lower  bound  on  the  error  probability  for  this  case  was  developed  in  [1], 
while  in  [2]  an  upper  bound  on  the  probability  of  bit  error  was  derived.  Phase  noise,  how¬ 
ever,  causes  each  linger  of  the  RAKE  receiver  to  have  the  same  tap  weight  estimation 
phase  error  at  any  instant  in  time,  since  the  received  pilot  tone  is  used  to  estimate  the  tap 
weight.  The  effects  of  phase  noise  on  RAKE  receivers  was  first  published  in  [3]  by  Zie- 
mer,  Vojcic,  Milstein,  and  Proakis;  however,  the  analysis  did  not  include  convolutional 
coding  or  pulse-noise  interference.  Both  commercial  and  military  systems  use  convolu¬ 
tional  coding  to  improve  the  bit  error  rate.  Military  systems  often  have  to  provide  reliable 
communication  in  the  presence  of  interference,  hence  it  is  important  to  extend  the  analy¬ 
sis  in  [3]  to  include  the  effect  of  coding  and  interference.  This  dissertation  will  follow  the 
analytical  approach  in  [3]  with  the  appropriate  modifications  for  pulse-noise  interference 
and  convolutional  coding.  In  our  work,  we  assume  that  a  pilot  tone  is  transmitted  on  a 
separate  frequency  and  is  an  integer  multiple  of  two  (2n)  times  the  carrier  frequency  of 
the  data  signal.  The  received  pilot  tone  then  serves  as  the  input  to  a  phase-lock-loop 
(PLL)  and  is  translated  down  to  the  carrier  frequency  of  the  data  so  it  can  be  used  for  co¬ 
herent  demodulation,  as  shown  in  Figure  1.3. 


Received  Pilot  Tone 


First  Order 

Divide  by 

Phase-Lock-Loop 

W 

2*n 

Pilot  Tone 


For  Demodulation 


Figure  1.3  Processing  Received  Pilot  Tone. 


The  analysis  is  done  for  various  levels  of  the  PLL  loop  SNR  in  order  to  determine  how 
coherent  systems  perform  for  different  values.  It  is  also  assumed  that  the  pilot  tone  is  re¬ 
ceived  with  fiat  fading  and  either  AWGN  or  AWGN  plus  a  noise-like  interference  signal 
that  is  modeled  as  additive  Gaussian  noise.  With  these  assumptions,  the  effect  of  pulse- 
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noise  interference  with  various  levels  of  received  signal-to-noise-ratios  on  the  pilot  tone 
can  be  investigated.  The  effect  of  AWGN  on  a  first  order  phase-locked  loop  is  a  Tikonov 
distributed  phase  error  as  shown  in  [4],  We  will  show  that  hard  decision  Viterbi  decoding 
can  perform  better  than  soft  decision  Viterbi  decoding  in  the  presence  of  Tikonov  distrib¬ 
uted  phase  noise. 

Commercially  available  soft  decision  Viterbi  decoders  are  designed  for  AWGN 
channels.  In  [5]  it  was  shown  that  these  decoders  perform  poorly  in  the  presence  of  pulse- 
noise  interference  and  that  hard  decision  Viterbi  decoding  performs  better  when  a  pulse- 
noise  interference  signal  is  present.  Several  modifications  to  improve  the  performance  of 
soft  decision  Viterbi  decoding  have  been  examined  in  the  literature.  In  [6]  the  concept  of 
a  noise-normalized  receiver  was  introduced,  and  simulation  results  in  [7]  show  that  the 
clipper  receiver  works  slightly  better  than  the  noise-normalized  receiver.  In  this  disserta¬ 
tion,  the  performance  of  a  maximum-likelihood  RAKE  receiver  optimized  for  both 
AWGN  and  pulse-noise  interference  is  derived,  and  the  resulting  performance  is  com¬ 
pared  with  those  obtained  using  the  RAKE  receiver  designed  for  AWGN  only. 

There  are  various  types  of  channel  fading,  but  in  this  dissertation  frequency- 
selective,  Ricean  distributed  fading  is  assumed.  Ricean  fading  has  two  components,  a  di¬ 
rect,  or  line-of-sight,  component  and  a  diffuse  component.  The  ratio  of  the  power  of  the 
direct  to  diffuse  component  will  be  referred  to  as  K.  When  A  =  0  there  is  no  direct 
signal  component,  and  this  is  referred  to  as  Rayleigh  fading.  If  A-^oo  then  there  is  no 
diffuse  component,  hence  no  fading,  which  happens  in  line-of-sight  communications 
when  there  is  nothing  between  the  transmitter  and  receiver  that  can  scatter  the  signal. 

Frequency-selective  fading  implies  that  different  spectral  components  of  the  sig¬ 
nal  are  affected  differently  by  the  channel.  For  flat  fading  channels,  all  the  frequency 
components  of  the  transmitted  signal  experience  the  same  amount  of  fading.  It  can  be 
shown  that  the  frequency-selective  channel  can  be  modeled  as  an  infinitely  long  tapped 
delay  line  with  tap  spacing  MW  where  W  is  the  noise-equivalent  bandwidth  of  the  base¬ 
band  signal,  and  time-varying  tap  weight  coefficients  [8].  For  a  given  channel  multipath 
spread,  we  know  that  for  time  delays  greater  than  the  multipath  spread  there  is  virtu¬ 
ally  no  signal  energy.  Hence,  a  very  good  model  for  the  frequency-selective  channel  is 
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obtained  by  using  a  tapped  delay  line  with  only  L  =  T^W  taps.  For  the  frequency- 
selective  channel  where  »  IjW ,  the  time  delay  IjW  between  successive  taps  is  suffi¬ 
cient  to  insure  that  each  of  the  tap  weights  are  mutually  uncorrelated  [8].  This  model  of  a 
frequency-selective  channel  is  shown  in  Figure  1.4. 


Figure  1.4  Model  of  the  frequency-selective  fading  channel. 


The  variables  «j ,  a^,  and  are  the  channel  tap  weights  that  scale  the  magnitude  of  the 

V"  signal  component  by  \a,\  and  the  phase  of  the  V’'  signal  component  by  Za;,and 

is  a  fixed  time  delay.  The  summer  adds  the  L  baseband  signals  together,  and  the  result  is 
a  frequency-selective  fading  signal. 

The  adjacent  bits  interfere  with  each  other,  as  shown  in  Figure  1.5,  when  several 
copies  of  the  transmitted  signal  are  received  with  different  delays.  If  the  three  signals 
shown  in  Figure  1.5  were  added  together,  the  interference  between  adjacent  bits  causes 
intersymbol  interference  (ISI),  which  increases  the  bit  error  rate.  In  a  coherent  RAKE  re¬ 
ceiver  the  L  different  signals  are  separated,  phase  shifted  to  remove  the  phase  offset  im¬ 
parted  by  the  communication  channel,  and  then  coherently  added  together. 
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Figure  1.5  Three  versions  of  the  same  signal  with  different  amplitudes  and  delays. 

This  eliminates  ISl  and  inereases  the  reeeived  signal-to-noise  ratio  (SNR).  Unfortunately, 
phase  noise  degrades  the  RAKE  receiver  channel  tap  weight  estimation  as  well  as  the  re¬ 
ceived  signal  demodulation.  In  [3]  Ziemer,  Vojcic,  Milstein,  and  Proakis  did  not  include 
both  of  these  effects. 

A  RAKE  receiver  separates  the  effective  L  different  received  signals  in  a  signal 
transmitted  over  a  frequency-selective  channel  by  using  the  correlation  properties  of  di¬ 
rect  sequence  spread  spectrum  (DSSS)  signals.  The  correlation  function  for  a  DSSS  sig¬ 
nal  is  1  for  zero  delay,  but  falls  linearly  to  -1/N ,  where  N  is  the  length  of  the  spreading 
sequence,  for  delays  of  one  chip  or  greater.  The  second  generation  digital  standard  for 
cellular  telephones  1S95  uses  a  sequence  that  has  a  length  of  2'^  -1  for  the 
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forward  channel  spreading.  As  a  result,  the  eorrelation  is  -  ij  «  0  for  delays  greater 

than  a  ehip.  This  means  that  any  multipath  eomponents  arriving  with  a  delay  greater  than 
one  ehip  will  have  zero  eorrelation. 

The  analysis  of  the  eoherent  RAKE  reeeiver  assumes  that  a  pilot  tone  and  the  data 
are  transmitted  on  separate  frequeneies  so  that  the  effeets  of  interferenee  on  the  earrier 
and  interferenee  on  the  data  ean  be  analyzed  separately.  The  pilot  tone  is  assumed  to  have 
a  fixed  average  SNR,  whieh  enables  the  random  phase  noise  effeets  to  be  analyzed  by 
fixing  the  random  phase  offset  and  then  integrating  over  the  probability  density  funetion 
(pdf)  of  the  phase  noise  to  remove  the  dependenee  on  the  random  phase  offset.  However, 
this  method  eannot  be  used  when  the  pilot  tone  phase  noise  is  time  varying,  espeeially 
when  the  pdf  for  the  phase  noise  depends  on  the  SNR  of  the  pilot  tone.  Commereial  fS95 
systems  transmit  the  pilot  tone  and  data  on  the  same  frequeney,  but  eneode  them  with 
separate  Walsh  funetions.  The  pilot  tone  is  eneoded  with  the  all  ones  Walsh  funetion, 
whieh  allows  the  reeeiver  to  reeover  the  pilot  tone  with  a  lowpass  filter,  fn  this  ease, 
pulse-noise  interferenee  would  be  able  to  affeet  the  pilot  tone  and  the  data  at  the  same 
time,  whieh  would  eause  a  time-varying  phase  noise  on  the  reeeived  pilot  tone,  and  it 
would  not  be  possible  to  derive  analytie  results.  Analytie  results  presented  here  are  ap- 
plieable  to  an  fS95-like  system  with  a  eonstant  phase  noise  level.  For  AWGN  ehannels 
without  phase  noise,  soft  deeision  Viterbi  deeoding  is  better  than  hard  deeision  Viterbi 
deeoding  [5];  however,  we  will  show  that  hard  deeision  Viterbi  deeoding  is  better  than 
soft  deeision  Viterbi  deeoding  when  the  pilot  tone  is  reeeived  with  a  10-dB  SNR  or 
lower. 

For  oases  of  severe  fading,  a  nonooherent  RAKE  reeeiver  is  generally  used  sinoe 
nonooherent  reoeivers  do  not  need  to  traok  the  phase  of  the  reeeived  signal.  This  means 
that  noncoherent  RAKE  reoeivers  are  not  affeoted  by  phase  noise  and  do  not  have  to  esti¬ 
mate  the  phase  of  the  ohannel  tap  weight.  The  drawbaok  of  using  a  nonooherent  RAKE 
reeeiver  is  that  the  reeeiver  will  have  nonooherent  oombining  losses  when 


7 


combining  noncoherently  detected  signals  [8].  The  performanee  of  a  noncoherent 
reeeiver  in  the  presenee  of  pulse-noise  interferenee  will  be  eompared  and  eontrasted  with 
that  of  a  eoherent  RAKE  receiver. 

The  coherent  and  noneoherent  RAKE  reeeivers  will  be  eompared  using  an  expo¬ 
nential  multipath  intensity  profile,  which  is  the  power  of  the  multipath  signal  as  a  func¬ 
tion  of  delay  for  a  given  instant  of  time.  The  performance  of  the  RAKE  receiver  is  best  for 
a  uniform  multipath  intensity  profile;  however,  real  world  systems  seldom  operate  over 
ehannels  with  uniform  multipath  intensity  profiles.  Sinee  an  exponential  multipath  inten¬ 
sity  profile  is  a  good  approximation  for  a  eongested  urban  area  [9],  in  this  dissertation  the 
multipath  intensity  profde  is  assumed  to  be  exp(-r^i^) ,  where  r  is  the  delay  and  y/  is 
determined  by  how  quickly  the  power  deereases  as  r  inereases.  In  this  work  we  assume 
that  the  first  finger  of  the  RAKE  receiver  has  a  delay  0  <  r  <  1 ,  the  seeond  finger  has  a  de¬ 
lay  1  <  r  <  2  the  third  finger  has  a  delay  2  <  r  <3,  and  so  on.  Both  the  coherent  and  the 
noncoherent  RAKE  receivers  will  be  analyzed  with  this  exponential  multipath  intensity 
profde. 

In  summary,  this  is  the  first  time  that  a  eoherent  RAKE  reeeiver  with  phase  noise, 
pulse-noise  interference  and  convolutional  coding  has  been  analyzed.  While  soft  decision 
Viterbi  decoding  outperforms  hard  decision  Viterbi  decoding  for  AWGN  channel  with 
perfect  phase  synchronization,  it  is  shown  that  the  reverse  is  true  when  the  pilot  tone  is 
received  with  a  10-dB  SNR  or  lower.  A  similar  result  was  obtained  in  [10]  where 
Shamain  and  Milstein  demonstrated  that  inereasing  the  antenna  diversity  beyond  two  ae- 
tually  increased  the  probability  of  bit  error  for  a  fading  ehannel  and  a  fixed  antenna  aper¬ 
ture  when  the  loop  SNR  was  5  dB.  In  order  to  get  the  probability  of  bit  error  to  eontinu- 
ally  deereased  as  the  antenna  diversity  increases,  Shamain  and  Milstein  found  that  a  loop 
SNR  of  20  dB  was  required.  In  [1 1,12]  performance  of  a  coded,  DSSS  signal  with  fading 
and  pulsed-noise  interference  was  analyzed,  but  these  results  did  not  include  AWGN  and 
phase  noise  effects.  The  maximum  likelihood  receiver  for  a  coherent  RAKE  receiver  with 
pulse-noise  interference  is  also  derived  in  this  dissertation  for  the  first  time.  Additionally, 
this  is  also  the  first  time  that  a  noncoherent  BFSK  RAKE  receiver  with  pulse-noise 
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interference  and  convolutional  coding  has  been  analyzed.  The  maximum-likelihood  re¬ 
ceiver  for  a  noncoherent  BFSK  RAKE  receiver  with  pulse-noise  interference  is  also  de¬ 
rived  and  analyzed. 

This  chapter  described  the  system  investigated,  presented  the  frequency-selective 
fading  channel  model  and  introduced  the  RAKE  receiver  and  imperfect  phase  synchroni¬ 
zation  concepts.  The  system  analyzed  was  compared  with  IS95-like  wireless  communica¬ 
tion  systems.  Finally,  related  research  was  discussed  with  an  emphasis  on  the  unique  con¬ 
tributions  from  this  dissertation. 

In  the  second  chapter,  the  coherent  RAKE  receiver  is  analyzed.  We  begin  by  ana¬ 
lyzing  the  coherent  RAKE  receiver  without  phase  noise  and  include  phase  noise  effects 
on  demodulation  and  channel-tap-weight  estimation  next.  Convolutional  coding  with 
both  hard  and  soft  decision  Viterbi  decoding  is  also  considered.  Ricean  fading,  pulse- 
noise  interference  and  AWGN  are  included  in  the  analysis.  In  order  to  optimize  the  per¬ 
formance,  the  maximum-likelihood  coherent  RAKE  receiver  for  pulse-noise  interference 
and  AWGN  is  derived  and  compared  with  the  performance  of  a  coherent  RAKE  receiver 
optimized  only  for  AWGN. 

Chapter  III  discusses  the  performance  of  noncoherent  RAKE  receivers  in  the  pres¬ 
ence  of  pulse-noise  interference  and  AWGN.  Once  again,  convolutional  coding  with  both 
hard  and  soft  decision  Viterbi  decoding  is  included  in  the  analysis.  The  received  signal  is 
assumed  to  experience  Ricean  fading,  pulse-noise  interference  and  AWGN.  The  maxi- 
mum-likelihood  receiver  for  pulse-noise  interference  and  AWGN  was  derived  for  the 
noncoherent  RAKE  receiver  and  the  performance  results  were  compared  with  a  nonco¬ 
herent  RAKE  receiver  designed  for  only  AWGN  only.  This  comparison  is  of  interest  be¬ 
cause  commercial  RAKE  receivers  are  designed  for  AWGN  only  channels. 

Next,  Chapter  IV  summarizes  the  findings  and  presents  recommendations  for  fur¬ 
ther  research. 


9 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


10 


II.  COHERENT  RAKE  RECEIVER  PERFORMANCE 

The  received  signal  obtained  after  transmission  over  a  frequency-selective  multi- 
path  channel  as  shown  in  Figure  1 .4  is  given  by 

r(t)  =  Re[f(t)e^"'],  (2.1) 

where  co  is  the  carrier  frequency  in  radians  and  the  complex  envelope  r{t^  is  given  by 

I 

?{t)  =  (2.2) 

;=1 

The  signal  r  (t)  can  be  recovered  by  a  RAKE  receiver  with  k  lingers,  as  shown  in  Figure 

2.1,  where  42 is  the  unfaded  carrier  amplitude,  a^  is  the  channel  tap  weight  deter¬ 
mined  by  the  multipath  intensity  profile, 

c(0  =  Re[c(t)e^'“'']  (2.3) 

is  the  locally  generated  chipping  waveform  where  c(t)  is  the  complex  envelope  of  the 


y{t)  =Re() 


Figure  2. 1  Coherent  RAKE  receiver. 
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chipping  waveform,  and  d{t)  is  the  baseband  data  waveform  which  equals  when  a 
binary  “1”  is  transmitted  and  when  a  binary  “0”  is  sent.  The  channel  is  assumed  to 
vary  slowly  eompared  to  the  bit  period  so  that  a.  can  be  modeled  as  a  constant  during 
each  bit.  Each  of  the  k  components  is  modeled  as  a  flat  fading  signal  and  can  have  a  dif¬ 
ferent  fading  coefficient  a^ ,  which  has  a  magnitude  and  a  phase  .  For  this  dis¬ 
sertation,  the  channel  tap  weights  are  modeled  as  Ricean  distributed  random  variables. 

The  ehipping  waveform  eomplex  envelope  c{t) ,  derived  from  the  ehipping  se¬ 
quence  ,  is  assumed  to  have  a  chip  rate  matched  to  the  RAKE  finger  spacing  .  The 
autocorrelation  function  for  the  ehipping  waveform  is  given  by  [13] 

fl  if  r  =  0 


.N 


i^0<T<N, 


(2.4) 


which  can  be  approximated  by 

=  (5'(r)  for -A^<  r  <  (2.5) 

and  is  periodic  with  a  period  of  N.  This  approximation  is  good  for  DSSS  signals  with 
very  long  PN  sequences  and  for  orthogonal  sequences.  The  finger  of  the  RAKE  re- 
eeiver  multiplies  the  received  signal  by  2y/2A^  |a^.  |  c  * (t  -  kT^ ) g  elimi¬ 

nates  the  phase  difference  between  the  fingers  so  that  each  finger  can  be  coherently 
summed  together  and  integrated  over  the  period  of  a  data  bit.  The  operation  of  multiply¬ 
ing  by  a  known  referenee  signal  and  integrating  is  referred  to  as  a  correlation  reeeiver 
[14].  Sinee  the  multipath  components  of  the  reeeived  signal  are  delayed  by  one  chip  or 
more,  from  (2.5)  we  see  that  only  the  multipath  component  with  zero  delay  eompared  to 
the  reference  signal  will  eorrelate  with  the  referenee  signal  and  the  other  multipath  eom- 
ponents  will  integrate  to  zero.  The  eomplex  envelope  of  the  signal  at  the  output  of  the 
RAKE  receiver  is  given  by 
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(2.6) 


t-T^  '«=1 

y.'^42A^\a\c{T -iT^)d{T -iT^)J^  ^"‘^'^dv, 

i-\ 

and  the  real  baseband  output  ean  be  found  by  taking  the  real  part  of  (2.6).  Eliminating  the 
terms  that  integrate  to  zero  in  (2.6),  we  get 

y{l)  =  2A;d{l)f^\a,\\  (2.7) 

7=1 

Henee,  the  RAKE  reeeiver  separates  the  multipath  eomponents,  removes  the  phase  offset, 
and  then  eoherently  eombines  them.  Sinee  the  reeeiver  eolleets  all  reeeived  signal  paths 
that  earry  the  same  information  on  separate  fingers,  its  aetion  is  similar  to  the  way  dirt 
eolleets  between  the  fingers  of  a  RAKE  when  one  RAKE?,  the  ground  in  a  garden  and  thus 
was  named  a  RAKE  reeeiver. 

A,  TAP  WEIGHT  ESTIMATION 

The  ehannel  tap  weights  a^,  a^,  . . .,  a,  seale  both  the  amplitude  and  phase  of  the 
reeeived  signal.  Sinee  the  ehannel  shifts  the  phase  by  Za, ,  the  RAKE  reeeiver  must  shift 
the  phase  on  eaeh  finger  by  -Za,  for  eoherent  eombining.  A  bloek  diagram  illustrating 

the  teehnique  for  estimating  the  ehannel  tap  weights  is  shown  in  Figure  2.2.  First,  the  re¬ 
eeived  signal  is  mixed  with  a  loeally  generated  ehipping  waveform  that  is  both  ehip  and 
sequenee  synehronized  with  the  reeeived  signal.  Mathematieally,  the  output  of  the  first 
mixer  on  the  k'^'  finger  ean  be  written 

(t)  =  2V24^ |«.| f  *  {t  - kT,  - T,)d{t - iT,  - T^)d{t  - iT^  -  ^2.8) 

i=\ 

The  above  result  is  delayed  by  one  bit  to  allow  for  demodulation  and  then  multiplied  by 
the  previous  bit  deeision.  At  the  output  of  the  seeond  mixer,  we  have 

y2(t)  =  2V24^|a.|f  *(t-k7; -T,)d{t -iT,-T,)d^  [t -kT, 

i=l 

The  lowpass  filter  then  removes  the  high  frequeney  terms,  and  from  (2.5)  the  output  of 
the  lowpass  filter  is  zero  when  i^k  .  Sinee  c*{t- kT^  - 7), ) c (t - kT^  - 7), )  =  1 ,  the  eom- 

plex  output  of  the  k'*  finger  is 
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y{t)  =  242 \a, \d4t- kT,  - T, ) .  (2.10) 

Since  BPSK  signaling  is  used,  d  (^)  is  either  +1  or  -1 ,  and  d^  [t)  always  +1 .  Thus, 
(2.10)  ean  be  simplified  to 

y(t)  =  2V24|a,|e^^“‘.  (2.11) 


Figure  2.2  Tap  Weight  Estimator  for  a  Coherent  RAKE  Reeeiver. 
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From  (2.1 1)  we  see  that  the  tap  weight  estimator  provides  estimates  of  both  the  magni¬ 
tude  of  the  ehannel  tap  weight  |a^|  as  well  as  the  phase  shift  Za^.  and  scales  these  esti¬ 
mates  by  ,  the  unfaded  carrier  amplitude. 


B,  PHASE  NOISE  MODEL 


Conventional  1S95  wireless  communication  systems  can  recover  the  pilot  tone  on 
the  forward  channel  by  using  a  low  pass  filter.  This  process  introduces  AWGN  which 
mixes  with  the  AWGN  on  the  data  signal  and  produces  multiplicative  noise  when  the  re¬ 
covered  pilot  tone  (plus  AWGN)  is  used  to  demodulate  the  data  bits.  In  order  to  remove 
the  noise  from  the  pilot  tone,  we  can  pass  the  received  pilot  tone  through  a  bank  of 
matched  filters  which  are  matched  to  different  frequencies.  The  frequency  of  the  pilot 
tone  can  then  be  estimated  by  interpolating  between  the  two  filters  with  the  highest  out¬ 
put.  This  previous  step  eliminates  noise  on  the  pilot  tone  but  can  produce  a  slight  fre¬ 
quency  offset  in  the  synthesized  pilot  tone,  and  frequency  offsets  in  the  pilot  tone  degrade 
the  data  demodulation.  Finally,  we  can  use  a  phase-locked  loop  to  track  the  pilot  tone; 
however,  signal  fading  and  AWGN  cause  phase  noise  at  the  phase-locked  loop  output. 

Phase  noise  is  the  random  phase  variation  of  the  carrier  and  causes  the  carrier  to 
deviate  from  a  perfect  cosine  wave.  It  is  shown  in  [4]  that  a  received  signal  and  AWGN 
applied  to  the  input  of  a  first  order  phase-locked  loop  causes  Tikonov  distributed  phase 
noise  which  is  a  function  of  the  loop  SNR.  The  pdf  for  the  Tikonov  distributed  phase  er¬ 
ror  is  [4] 


exp(y0  cos(^zi)) 


(2.12) 


where  (5  is  the  loop  SNR,  which  can  be  affected  by  channel  fading,  and  1^  is  the  modi¬ 
fied  Bessel  function  of  the  first  kind  and  order  zero.  A  plot  of  the  Tikonov  pdf  for  various 
values  of  loop  SNR  is  shown  in  Figure  2.3.  The  solid  line  shows  the  Tikonov  pdf 
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(|)  =  -7t  to  7t 

Figure  2.3  F’(^|y0)  vs.  (j)  for  Tikonov  phase  noise. 

for  a  loop  SNR  of  1 5  dB,  and  is  similar  to  a  Gaussian  pdf.  The  trace  with  large  dashes  is 
for  a  loop  SNR  of  10  dB,  which  also  looks  like  a  Gaussian  pdf  The  trace  with  small 
dashes  is  for  a  loop  SNR  of  5  dB,  and  the  probability  of  large  phase  offsets  increases  for 
this  level  of  loop  SNR.  Recall  that  there  is  no  signal  present  when  the  loop  SNR  is  zero 
(-00  dB) ,  and  we  expect  the  phase  error  from  the  phase-locked  loop  to  be  uniformly  dis¬ 
tributed  from  -7t  to  TV  as  shown  in  Figure  2.3. 

C.  EFFECT  OF  PHASE  AND  FREQUENCY  OFFSETS  ON  BPSK  SIGNALS 

The  next  step  is  to  analyze  how  phase  noise  affects  the  received  signal.  Figure  2.4 
shows  a  BPSK  receiver  with  an  unknown  frequency  offset  Ao  and  phase  offset  (j)  in  the 
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local  oscillator  where  the  output  P{T)  is  assumed  to  be  sampled  every  T  seeonds.  The 
output  ean  then  be  written  as: 

2,\f2A  ^ 

P{T)  = - ^  jsm((»ot)sm[((yo  +  +  (f^dt.  (2.13) 

2"  0 


yflA^  sin(n>nt) 


2sm[((U(,  +  A(»)t  +  ^zJ] 

Figure  2.4  BPSK  reeeiver  with  unknown  frequeney  and  phase  offset. 
Substituting  the  identity 


into  (2.13),  we  obtain 


sm(a)  sm(h)  =  ^  [cos(a  -b)-  cos{a  +  h)] 


yflA 


(2.14) 


=  — ^|cos[(A(y)t  +  ^]  -eos[(2(»o  +  A(»)t  +  ^/i]  Jt.  (2.15) 

2"  0 

Equation  (2.15)  has  a  high  frequency  term  at  twice  the  carrier  frequency  and  a  baseband 
component.  Since  every  complete  cycle  of  the  high  frequency  term  will  integrate  to  zero, 
the  high  frequency  component  will  either  integrate  to  zero  or  a  very  small  number  if  there 
is  not  an  integral  number  of  eycles  over  the  bit  period.  Fienoe,  the  double  frequency  term 
is  negligible  and  only  the  baseband  term  will  be  present  at  the  output  of  the  integrator. 
Using  the  identity 

eos(a  +  h)  =  cos(a)eos(h)-sm(a)sm(h)  (2.16) 


in  (2.15),  we  get 

P{T)  = 


42A, 


1  1 

j  cos[(A(y)t]cos(^zi)it- jsm[(A(»)t]sm(^/J) 
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(2.17) 


which  can  be  integrated  to  produce 

P{T)  =  V24  cos(^)  r cos(A(ur)-l1.  (2.18) 

At  this  point  we  will  limit  our  diseussion  by  only  eonsidering  the  phase-locked  loop 
tracking  problem,  which  has  no  frequency  offset  but  unwanted  phase  noise  due  to  the  pi¬ 
lot  tone  being  received  with  AWGN.  Setting  A®  to  zero  in  (2.18),  we  get  the  result 

P(r)  =  V24cos(^/i).  (2.19) 

Hence,  the  effect  of  an  unknown  phase  error  at  the  BPSK  receiver  is  to  reduce  the  ampli¬ 
tude  of  the  carrier  by  the  cosine  of  the  phase  error. 

D,  EFFECT  OF  PHASE  NOISE  ON  CHANNEL  TAP  WEIGHT 
ESTIMATION 

Not  only  does  phase  noise  affeet  the  data  demodulation  in  a  RAKE  receiver,  it 
also  degrades  the  RAKE  receiver  channel  tap  weight  estimation.  We  begin  with 
the  tap  weight  estimator  shown  in  Figure  2.2  to  include  a  random  phase  offset  (j)  as 
shown  in  Figure  2.5.  For  a  signal  transmitted  over  a  frequeney-seleetive,  multipath  ehan- 
nel,  the  complex  envelope  of  the  received  signal  can  be  written  as  shown  in  (2.2).  In  this 
case,  however,  the  local  oscillator  in  the  receiver  is  slightly  out  of  phase  with  the  received 
signal,  so  the  locally  generated  chipping  sequence  includes  a  random  carrier  phase  offset 
(j) .  We  ean  write  the  eomplex  envelope  of  the  locally  generated  ehipping  sequenee  on  the 

k‘^  linger  as 

/  {t)  =  c*{t- kT, (2.20) 
Multiplying  (2.20)  by  the  complex  envelope  of  the  received  signal  (2.2),  we  get 

y,{t)  =  2V24  Y.  h  )c{t  -  iT,  )d{t  -  iT,  y-Ko^T,-za,)^-j(cokT,-^)  _  ^2.21) 

i=l 

The  above  result  is  delayed  by  one  bit  to  allow  for  demodulation  and  then  multiplied  by 
the  previous  bit  deeision.  At  the  output  of  the  seeond  mixer,  we  have 

72  (0  =  2V24  Y  I  c  -  ^2;  )d{t  -  iT,  )d{t-kT,)d{t-  iT,  _  ^222) 
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The  lowpass  filter  then  removes  the  high  frequency  terms,  and  from  (2.5)  the  output  of 
the  lowpass  filter  is  zero  when  i^k.  Again,  since  that  c*{t -kT^^c{t -kT^^  =  \ ,  the 

complex  output  of  the  k‘^  finger  is 


Figure  2.5  Channel  Tap  Weight  Estimator  with  Phase  Noise  (j) . 


y{t)  =  2V24  \a,\d\t-  kT, )  (2.23) 

Since  BPSK  signaling  is  used,  d  (t)  is  either  +1  or  -1 ,  and  d^  (t)  is  always  +1 .  Thus, 
(2.23)  can  be  simplified  to 

y(t)  =  2V24|a,|e^'(^“‘"^),  (2.24) 
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which  differs  from  the  tap  weight  estimate  without  phase  noise  only  by  the  term.  It 
is  interesting  to  note  that  phase  noise  affects  both  the  tap  weight  estimate  and  the  signal 
demodulation  and,  to  this  author’s  knowledge,  no  other  researeh  has  examined  the  effeets 
of  both  simultaneously. 

E,  PULSE-NOISE  INTERFERENCE 

In  addition  to  phase  noise,  interferenee  is  also  eonsidered.  The  interferenee  signal 
is  modeled  as  AWGN.  Two  separate  interference  scenarios  are  analyzed  since  the  trans¬ 
mitter  is  assumed  to  send  a  pilot  tone  on  a  separate  frequeney  than  the  data.  First,  inter¬ 
ference  of  only  the  data  signal  frequency  is  analyzed.  Seeond,  interference  of  the  pilot 
tone  is  analyzed.  When  the  pilot  tone  experiences  interference,  all  the  data  bits  are  af¬ 
fected  by  inereased  phase  noise  due  to  the  interference  since  the  pilot  tone  is  used  to  syn- 
ehronize  the  receiver  local  oscillator  in  order  to  demodulate  the  data  signal. 

Pulsing  the  interferenee  signal  is  also  investigated.  Speeifieally,  if  the  interference 
power  speetral  density  is  Nj  and  the  duty  eyele  is  p ,  then  the  interferenee  power  is 

Njl pTj,  where  is  the  bit  duration.  This  assures  a  eonstant  average  power  interferenee 
signal  regardless  of  p .  Pulsing  can  be  an  effeetive  interference  strategy  beeause  of  the 
inereased  power  when  interferenee  is  present.  Sinee  the  probability  of  bit  error  is  given 
by  =  /?  (error  rate  with  pulse-noise  interference)  +  (l  -  p)  (error  rate  with  AWGN), 
assuming  the  probability  of  bit  error  is  small  when  only  AWGN  is  present,  the  overall 
probability  of  error  is  approximately  5  x  10  ^^  if  the  interference  can  cause  a  50  percent 
probability  of  error  when  it  is  on  and  if  the  interferenee  duty  eyele  ( /?)  is  one  pereent. 

F.  MAXIMUM-LIKELIHOOD  RECEIVER  FOR  CHANNELS  WITH 

GAUSSIAN  NOISE,  PULSE  INTERFERENCE,  AND  DIVERSITY 

The  maximum-likelihood  receiver  for  a  BPSK  signal  with  d  independent  diver¬ 
sity  receptions,  of  whieh  i  have  pulse-noise  interference,  is  derived  in  this  seetion.  The 
reeeived  signal  is  designated  (t)  where  m  takes  values  between  1  and  d  to  represent 
eaeh  of  the  d  different  diversity  reeeptions.  The  interferenee  waveform  is  modeled  as 
additive  Gaussian  noise  with  varianee  <jI  when  only  AWGN  is  present  and  with  vari- 
anee  <jI  +  cr^  when  interferenee  is  present.  We  begin  with  a  general  analysis  that  is 
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valid  for  diversity  in  space,  frequency  and  time.  Assuming  that  the  diversity  receptions 
are  independent,  we  obtain  the  joint  pdf  for  the  “0”  bit  when  i  diversity  receptions  ex¬ 
perience  interference  and  d  -i  diversity  receptions  have  only  AWGN  as 


d-i 


Inal  j 


exp 


-1 


d-i 


2a, 


0  fn=l 


(2.25) 


1 


2;z-(cro  +a^)  j 


exp 


-1 


2(0-0  +^,) 


=d-i+\ 


dr, 


and  the  joint  pdf  for  the  “1”  bit  is  given  by 


fAy)=\ 


d-i 

1 


1 


exp 


2a, 


0  m=\ 


(2.26) 


2n{aa  +CTj)j 


exp 


-1 


2(0-0 


0  '  ^  j  )  m=d-i-\-\ 


dr. 


Forming  the  likelihood  ratio  [15],  we  obtain 

ZM>i 

/o  (k) 


(2.27) 


If  the  likelihood  ratio  is  greater  than  1,  the  receiver  decides  bit  “1”  was  sent,  and,  if  the 
likelihood  ratio  is  less  than  1,  the  receiver  decides  bit  “0”  was  sent.  Substituting  (2.25) 
and  (2.26)  into  (2.27)  and  simplifying,  we  get 


J  exp  ^^(2y„^(r)5o(r)-2y„,(r)5i(r)-5o(r)  +  5f(r)) 


t-T 


xexp 


0  »t=l 

-1 


Z  (2k.(^)A(^)-2k„(r)5i(r)-5o(r)  +  5f(r)) 


m=d-i-\-\ 


dr^l. 


(2.28) 


< 

“50 


Taking  the  natural  log  of  both  sides  and  rearranging  terms,  we  obtain 

-1  ^ 


^Z  1  {ymi^)Soi^)-y„,{T^)s,{T^))dT 


+ 


cr, 


0  m=\ 


2  I  2 

^0  '  ^ j  m=d-i+\ 


<  2a^ 

So  ^^0  »I=1  t-T 


+ 


-1 


2( 


2  2  \ 

^0  ^ j  )  i+1  f-X 


Z  1 


(2.29) 
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For  BPSK  signals,  note  that  (t)  =  -5j  (t)  so 

^0  (0-^1  (0  =  0-  (2.30) 

Substituting  (2.30)  into  (2.29),  we  get  the  final  form  for  the  log  likelihood  ratio  for  a  co¬ 
herent  RAKE  receiver  with  AWGN  and  pulse-noise  interference: 


Z  J  2]  2  Z  1  (k;„(^K(^)-k;«(^>i(^))^o.(2.3i) 

m=l  i-E  ^0  '  ^ j  fn=d-i+\  f-j  Sq 


An  implementation  of  the  maximum-likelihood  receiver  is  shown  in  Figure  2.6,  where 
l/cr^  represents  the  inverse  of  the  received  signal  variance,  which  is  equal  to  l/ cr^  when 

AWGN  is  present  and  l/(cro  +cr^) 


<j 


Figure  2.6  Maximum-Likelihood  Receiver  for  Diversity  Receptions  of  Baseband 

BPSK  Signals  with  Pulse  Interference. 

when  the  pulse-noise  interference  is  present.  The  receiver  structure  shown  in  Figure  2.6 
can  be  simplified  for  BPSK  signals.  Since  =  -5j ,  the  result  from  correlating  with  5g  is 

the  negative  of  the  result  from  correlating  with  ,  which  means  that  only  one  correlator 
is  necessary.  A  diagram  of  the  maximum  likelihood  RAKE  receiver  for  AWGN  and 
pulse-noise  interference  is  shown  in  Figure  2.7,  where  l/cr^  represents  the  inverse  of  the 

variance,  equal  to  l/cTg  when  only  AWGN  is  present  and  l/(cro  )  when  the  pulse- 
noise  interference  is  present. 
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y{t)  =  Re() 


Figure  2.7  Maximum-Likelihood  Coherent  RAKE  Receiver  for  AWGN  and  Pulse- 

Noise  Interference. 

This  RAKE  receiver  is  similar  to  the  one  shown  in  Figure  2.1,  except  that  the  out¬ 
put  is  multiplied  by  the  inverse  of  the  total  variance,  not  just  the  inverse  of  the  variance 
due  to  AWGN.  The  multiplication  by  the  inverse  of  the  variance  complicates  the  RAKE 
receiver  because  the  variance  must  be  computed  for  every  bit,  which  also  requires  addi¬ 
tional  circuitry  to  compute  the  variance. 

G.  MEAN  AND  VARIANCE  OF  COHERENT  RAKE  RECEIVERS 

We  stated  earlier  that  pulse-noise  interference  causes  a  time-varying  variance 
which  we  now  investigate  further.  Additionally,  we  present  the  method  of  analysis  for 
coherent  R/IAE’  receivers  with  pulse-noise  interference  and  convolutional  coding.  First, 
we  begin  by  assuming  that  the  received  signal  can  be  modeled  as  a  Gaussian  random 
variable  with  a  conditional  mean  of  yflA^  |«^.| ,  where  |a^|  is  the  random  Ricean  faded  tap 
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weight,  and  a  variance  of  .  Since  |a^  |  is  modeled  as  a  random  variable,  we  obtain  the 
conditional  probability  of  bit  error  and  then  average  over  the  pdf  of  the  random  variable 
|a^.|  to  get  an  average  probability  of  bit  error.  The  first  mixer  in  Figure  2.1  multiplies  the 

received  signal  by  the  chipping  sequence,  which  removes  the  spreading  but  does  not 
change  the  probability  distribution  of  the  received  signal.  The  second  mixer  in  Figure  2.1 

multiplies  the  down-converted  received  signal  by  This  scales  the  signal 

by  V24 1  «^|  and  phase  shifts  the  multipath  components  by  for  coherent  combining. 
Hence,  the  output  from  the  second  mixer  in  Figure  2.1  is  a  Gaussian  random  variable 
with  mean  2A^  and  variance  2cr  A^  .  The  coherent  receiver  in  Figure 
2.1  then  sums  together  the  j  fingers.  This  is  the  point  where  pulse  interference  compli¬ 
cates  the  analysis.  If  only  AWGN  is  present,  then  the  mean  at  the  output  of  the  coherent 
RAKE  receiver  can  be  written  as 

m  =  2AlY\^a,\\  (2.32) 

and  the  variance  at  the  output  of  the  RAKE  receiver  is 

al=2a-Aljy,,\\  (2.33) 

./=1 

where  cr^  =  cr^  when  only  AWGN  is  present  and  cr^  =  erg  +  crj  when  the  pulse-noise  in¬ 
terference  is  present.  We  require 

ZKf=l.  (2.34) 

k=\ 

The  BPSK  bit  error  rate  probability  for  the  RAKE  receiver  is  given  by 

( 

—  =Pr(y>0|0)Pr(0)  +  Pr(y<0|l)Pr(l),  (2.35) 

\^r) 

where  y  is  the  random  variable  that  represents  the  RAKE  output.  Due  to  the  symmetry  of 
both  the  receiver  and  the  noise,  wheny  is  modeled  as  a  conditional  Gaussian  random 
variable,  (2.35)  can  be  evaluated  to  obtain  the  conditional  probability  of  bit  error  [17] 
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Substituting  (2.32)  and  (2.33)  into  (2.36),  we  get 


Equation  (2.37)  ean  be  simplified  to  obtain 


By  defining 


yb  =  — , 
a 


(2.36) 


(2.37) 


(2.38) 


(2.39) 


the  probability  of  bit  error  eonditioned  on  ean  be  written 

(2.40) 

This  is  equation  (14.3-2)  on  page  817  of  [9],  whieh  shows  how  to  remove  the  depend- 
enee  on  for  the  ease  of  Rayleigh  distributed  fading. 

Pulse  interferenee  eomplieates  the  analysis  sinee  the  interferenee  ean  be  on  for 
only  a  fraetion  of  a  bit,  in  whieh  ease  eaeh  finger  would  have  a  different  varianee.  In  this 
ease  (2.33)  has  to  be  written  as 

(2.41) 

k=l 

where  al  is  cTq  when  only  Gaussian  noise  is  present  and  erg  +  <j^  when  the  pulse-noise 

interferenee  is  present.  To  prevent  this  problem,  we  assume  that  eaeh  bit  will  either  ex- 
perienee  interferenee  for  the  full  bit  period  or  no  interferenee.  Sinee  we  are  not  eonsider- 
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ing  cases  where  a  fraetion  of  a  bit  experienees  interferenee,  (2.36)  is  valid,  with  the  addi¬ 
tional  eonsideration  that  is  eonditional  on  whether  or  not  the  bit  experienees  interfer¬ 
enee. 

The  analytie  approaeh  taken  here  is  to  obtain  the  pdf  at  the  output  of  eaeh  finger 
of  the  RAKE  reeeiver  and  eonvolve  them  together  to  get  the  pdf  for  the  random  variable 
that  represents  the  sum  of  the  RAKE  finger  outputs.  We  perform  the  analysis  assuming 
that  a  bit  “0”  was  transmitted.  The  probability  of  bit  error  is  then  determined  by  integrat¬ 
ing  the  pdf  of  the  random  variable  that  represents  the  output  of  the  RAKE  reeeiver  from  0 
to  00 .  The  eonvolutions  are  performed  by  taking  the  two-sided  Laplaee  transform  of  eaeh 
finger’s  pdf,  multiplying  the  resulting  Laplaee  transforms  together,  and  then  taking  the 
inverse  Laplaee  transform,  whieh  when  integrated  from  0  to  oo  gives  the  probability  of 
bit  error. 

H.  ANALYSIS  OF  A  COHERENT  RAKE  RECEIVER  WITH  PHASE  NOISE 
AND  RICEAN  FADING 

In  Figure  2.1  the  eoherent  RAKE  reeeiver  shown  does  not  inelude  phase  noise. 
We  have  already  seen  that  phase  noise  eauses  a  random  phase  error  (j)  between  the  lo- 
eally  generated  ehipping  sequenee  and  the  reeeived  signal  and  that  the  eomplex  baseband 
signal  is  multiplied  by  .  The  reeeived  signal  shown  in  Figure  2.8  now  ineludes  the 
AWGN  n(t)  =  where: 

I 

r(t)  =  n(t)  +  ^  yflA^afiit  -  iT,  )d{t  -  iT,  .  (2.42) 

!=0 

Reeall  that  multiplying  the  reeeived  signal  by  a  ehipping  sequenee  does  not  affeet  the 
Gaussian  noise,  so  at  the  output  of  the  first  mixer  we  get 

y,{t)  =  h{t)  +  ly/lA^  ^  |«.  I  f  *  (t  -  kT, )d{t  -  iT, )d{t  -  iT, -*) .  (2.43) 

;=0 

whieh  then  is  multiplied  by  the  ehannel  tap  weight  -JlA^  ,  summed  with  the 

output  of  the  other  fingers,  and  integrated.  The  real  baseband  output  is  then  found  by  tak¬ 
ing  the  real  part  of  the  integrator  output 

y{t)  =  2Al  cos  (2^)  J  (t)  ^  |a .  |  V  n  (t) .  (2.44) 

y=i 
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y{t)  =Re() 
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Equation  (2.44)  is  the  output  for  the  summation  of  all  the  RAKE  linger  outputs.  The  mean 
value  at  the  output  of  the  finger  is 

m  =  2 A]  cos  (2^)  ,  (2.47) 

and  the  variance  is 

a]  =2A^^zo%{2(f)y^(7^ ,  (2.48) 

where  =  crl  when  only  AWGN  is  present  and  cr^  =  (Tq  +  cr]  when  the  pulse  interfer¬ 
ence  is  present.  From  (2.47)  and  (2.48)  we  see  that  both  the  mean  and  the  variance  at  the 
E’"  finger  output  are  dependent  on  the  noncentral  chi-squared  distributed  as  well  as 
the  Tikonov  distributed  (j) .  Substituting  (2.47)  and  (2.48)  into  (2.46),  we  get 


-24V,  cos(2^/i) 


y 

2_2  A 

s  a 

s  - 

V 

2  JJ 

(2.49) 


We  now  remove  the  dependence  on  by  integrating  over  the  chi-squared  distribution 
function.  We  rewrite  (2.49)  as 


where 

p  =  2Alcos{2(l>) 

Using  (2.45)  and  (2.50),  we  obtain 

00  Y 

(■^  I  ^^^)  =  f  exp  )  tr^exp 

n  ^O'k 


s - 


V 


(2.50) 


(2.51) 


J 


-\  +  K 

J 


cr, 


k  ’ 


(2.52) 


k  J 


which  can  be  written 

00  -y 

O,  (5  I  =  exp(-A)J  — exp 


2at 


2Ky, 


kJA  V 


cr 


dy,. 


(2.53) 


k  J 


Equation  (2.53)  can  be  evaluated  using  the  identity  (6.643.4)  from  [18] 


(2/34^yx  =  n ! /3''e  “ 


V  «  y 


(2.54) 
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where  (x)  is  the  Lagerre  polynomial.  The  Bessel  funetion  of  the  first  kind  and  order 
zero  Jq  is  related  to  the  modified  Bessel  funetion  of  the  first  kind  and  order  zero  Iq  by 
/q (x)  =  Jq{ jx)  [18]  where  j  =  .  From  (2.53)  we  see  that  n  =  v  =  0 ,  henee  we  ean 

use  the  relationship  [18] 

T:(x)  =  1.  (2.55) 


Therefore,  we  obtain  the  Laplaee  transform  of  the  pdf  at  the  output  of  the  k‘^  finger  as 

1 


0,(5|^)  = 


lalp  +  l 


exp 


^  -IK 


(2.56) 


where  fi  is  defined  in  (2.51)  and  2cr^  is  the  diffuse  power  in  eaeh  finger.  The  RAKE  re- 
eeiver  adds  the  outputs  of  the  fingers  together,  so  the  Laplaee  transform  of  the  pdf  of  the 
random  variable  that  represents  the  output  of  the  RAKE  reeeiver  ean  be  written 


(2.57) 


k=\ 


To  remove  the  dependence  on  ^  we  must  still  average  over  the  pdf  of  (/) : 

=  (2.58) 

Equation  (2.58)  must  be  evaluated  numerically. 

The  probability  of  bit  error  is  the  integral  from  0  to  oo  of  the  inverse  Laplace 
transform  of  d)  (5) .  The  inverse  Laplace  transform  and  the  subsequent  integration  can  be 

done  recursively  using  the  Gauss-Chebyshev  quadrature  method.  Defining  Y  as  the  ran¬ 
dom  variable  at  the  output  of  the  RAKE  receiver  with  a  pdf  given  by  the  inverse  Laplace 
transform  of  d)(^)  we  have  [3] 


1  nil 

P{Y  >0)  =  — ^Re[o(5  =  c  +  ycr^.)]  +  Im[o(5  =  c  +  ycr^)]  -t-L,  (2.59) 

n  k=\ 


where 

=  tan[(2A:-l);r/2n],  (2.60) 

c  is  a  constant  chosen  to  make  the  algorithm  converge  quickly  and  E  is  an  error  term 
that  approaches  0  as  n  gets  large.  Extremely  low  bit  error  rates  require  larger  values  of 
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n  .  Values  of  n  between  200  and  800  were  used.  Methods  for  ehoosing  the  parameter  c 
are  diseussed  in  [19],  and  examples  of  MATLAB  eode  to  implement  (2.59)  are  given  in 
Appendix  A. 

1,  Addition  of  Pulse  Interference 

The  preeeding  analysis  does  not  inelude  the  effeet  of  pulse-noise  interferenee,  but 
it  ean  be  easily  modified  to  do  so.  It  is  assumed  that  the  varianee  of  the  reeeived  noise  is 
cTq  when  only  AWGN  is  present  and  cr^  +  cr^  when  the  pulse  interferenee  is  present.  The 

only  term  in  the  preeeding  analysis  that  depends  on  the  varianee  of  the  reeeived  signal  is 
We  now  define 


(5'j  =  2AI  cos{2(/)) 


S  <7, 


S-- 


-2  A 
0 


(2.61) 


and 


5 2  =  2AI  cos{2(/)) 


s - 


(2.62) 


V  J 

When  the  reeeived  signal  has  only  AWGN,  we  simply  substitute  f3  =  5^  into  (2.56)  and 


when  the  interferenee  signal  is  present  we  substitute  ^  =  3^  into  (2.56),  leading  to 

f  I  ^ 


V  1=1 


(2.63) 


J 


when  only  AWGN  is  present  and 


0(5|^)  = 


(2.64) 


V  1=1  J 

when  all  the  fingers  experienee  interferenee.  It  should  be  mentioned  that  (2.63)  and  (2.64) 

are  conditional  on  (j)  and  must  be  numerically  integrated  over  the  Tikonov  pdf  to  remove 

the  conditioning  on  (j)  as  shown  in  (2.58). 

2.  Maximum-Likelihood  RAKE  Receiver  for  AWGN  and  Pulse-Noise  In¬ 
terference 

The  maximum-likelihood  RAKE  receiver  for  AWGN  and  pulse-noise  interference 
shown  in  Figure  2.7  has  an  additional  mixer  that  multiplies  the  received  signal  by  the  in- 
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verse  of  the  variance.  Since  it  is  assumed  that  the  bit  either  completely  experiences  inter¬ 
ference  or  no  interference,  this  receiver  multiplies  the  received  signal  by  l/ when  only 

AWGN  is  present  and  by  ^  when  the  pulse-noise  interference  is  present.  Mul¬ 

tiplying  by  the  inverse  of  the  variance  affects  both  the  mean  and  the  variance  of  the  signal 
at  the  output  of  the  RAKE  linger.  In  this  case,  the  mean  is  given  by  (2.47)  divided  by  <J^ 

_  2A^ 

(2-65) 

cr 

and  the  variance  is  given  by  (2.48)  divided  by 

(2-66) 


where  cr^  is  gI  when  only  Gaussian  noise  is  present  and  gI  +  g^-  when  the  pulse-noise 

interference  is  present.  The  variance  of  a  RAKE  receiver  designed  for  AWGN  channels  is 
proportional  to  the  variance  of  the  received  signal  as  shown  in  (2.48).  For  the  maximum- 
likelihood  RAKE  receiver  designed  for  AWGN  and  pulse-noise  interference,  both  the 
mean  and  variance  are  proportional  to  the  inverse  of  the  variance.  This  means  that  a  high 
powered  interference  signal  will  tend  to  drive  the  RAKE  receiver  output  toward  zero,  as 
opposed  to  the  RAKE  receiver  designed  only  for  AWGN  where  a  high  powered  interfer¬ 
ence  signal  will  drive  the  output  towards  infinity.  The  analysis  for  this  maximum- 
likelihood  RAKE  receiver  is  the  same  as  the  previous  analysis  except  that  in  this  case 


= — Acos(2^zi) 


CTn 


when  only  AWGN  is  present  and 

.  24' 

^4  = 


cos 


(2^) 


V 


1- 


(2.67) 
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■^(^0 


(2.68) 


when  the  interference  signal  is  present.  Once  again  we  get 


f  I 


when  only  AWGN  is  present  and 


V  !=1 


(2.69) 
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(2.70) 


^[s\(l),p  =  d,)=  Y\^,[s\(l),p  =  d,) 

V  1=1  y 

when  all  the  lingers  experienee  interferenee.  The  result  is  numerieally  integrated  over  the 
Tikonov  pdf  to  remove  the  dependenee  on  (j)  as  shown  in  (2.58). 

3,  Hard  Decision  Decoding 

The  previous  analysis  eomputed  the  uneoded  bit  error  rate.  Commereial  wireless 
standards  sueh  as  IS95  use  a  eonvolutional  eode  to  improve  the  bit  error  ratio.  In  this  see- 
tion  we  analyze  the  performanee  of  a  eonvolutionally  eneoded  signal  transmitted  over  a 
Rieean  distributed,  frequeney-seleetive  fading  ehannel  with  pulse-noise  interferenee  and 
hard  deeision  Viterbi  deeoding. 

The  probability  of  bit  error  is  upper  bounded  by  [9] 

00 

(2.71) 

‘l=‘l  free 

where  is  the  total  number  of  nonzero  information  bits  in  all  weight  d  paths,  is  the 
probability  of  seleeting  a  eode  sequenee  that  is  a  Hamming  distanee  d  away  from  the  eor- 
reet  eode  sequenee,  and  d^^^^  is  the  minimum  weight  path  in  the  eonvolutional  code  state 

diagram  that  diverges  from  and  remerges  with  the  all-zero  state  [20].  For  hard  decision 
Viterbi  decoding,  P^  is  given  by  [9] 

2 


(2.73) 


(2.74) 
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over  the  Tikonov  pdf,  numerically  invert  the  Laplace  transform,  and  integrate  the  result 
from  0  to  00 .  In  (2.74)  p  is  the  duty  cycle  of  the  interference,  5^  is  defined  in  (2.61)  and 
is  defined  in  (2.62).  Figure  2.9  shows  the  probability  of  bit  error  for  a  three  finger 
RAKE  receiver  with  rate  14,  constraint  length  9  convolutional  coding  and  hard  hard  deci¬ 
sion  Viterbi  decoding  where  =  15  dB,  the  pilot  tone  has  a  20-dB  loop  SNR,  the 

ratio  of  direct-to-diffuse  power  is  10  dB,  and  the  interference  signal  has  a  one  percent 
duty  cycle  ( /?  =  0.01 ),  a  10%  duty  cycle  ( /?  =  0. 1 ),  and  a  100%  duty  cycle  (/?  =  !).  For  a 
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Figure  2.9  Coherent  RAKE  Receiver  with  Flard  Decision  Viterbi  Decoding. 

rate  14,  constraint  length  9  convolutional  code  the  first  nonzero  three  terms  of  the  series  in 
(2.71)  are  337^2 +2817^4 +2179/^g .  An  exponential  multipath  intensity  profile  (MIP)  is 

assumed  where  each  finger  has  10  dB  less  signal  energy  than  the  previous  finger.  We  see 
from  Figure  2.9,  that  the  interference  should  be  on  all  the  time  if  it  has  approximately  the 
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same  power  as  the  reeeived  signal  when  hard  deeision  Viterbi  deeoding  is  used.  When 
the  interferenee  has  less  power  than  the  reeeived  signal,  then  pulsing  the  interferenee  sig¬ 
nal  is  the  best  strategy.  It  is  interesting  to  note  that  the  error  rate  is  extremely  low  for  all 
values  of  E^jNj  when  the  interferenee  duty  cyele  is  0.01.  We  will  show  that  pulsing  the 

interferenee  signal  is  mueh  more  effeetive  against  soft  deeision  Viterbi  deeoding. 

For  the  ease  of  the  maximum-likelihood  RAKE  reeeiver  designed  for  AWGN  and 
pulse-noise  interferenee,  the  probability  of  making  an  error  p  ean  be  ealeulated  by  inte¬ 
grating 


0(5|^)  =  p(0(5|^/i,J  =  J3))  +  (l-p)(0(5|^,J  =  J4))  (2.75) 

over  the  Tikonov  pdf  as  shown  in  (2.58),  numerieally  inverting  the  Laplaee  transform  and 
integrating  the  result  from  0  to  go  .  The  results  for  hard  deeision  Viterbi  deeoding  using 
the  maximum-likelihood  RAKE  reeeiver  designed  for  AWGN  and  pulse-noise  interfer¬ 
enee  were  the  same  as  the  RAKE  reeeiver  designed  for  AWGN.  This  ean  be  explained  by 
noting  that,  for  the  RAKE  reeeiver  designed  for  AWGN,  the  interferenee  signal  eauses  the 
varianee  of  the  reeeived  signal  to  be  very  large  when  the  interferenee  signal  is  on,  whieh 
in  turn  ean  lead  to  a  50  pereent  error  rate.  The  maximum-likelihood  RAKE  reeeiver  de¬ 
signed  for  AWGN  and  pulse-noise  interferenee  eauses  the  mean  and  varianee  to  tend  to¬ 
wards  zero  when  the  interferenee  signal  is  present.  Sinee  hard  deeision  Viterbi  deeoding 
makes  deeisions  on  a  bit-by-bit  basis,  a  reeeived  signal  with  a  mean  and  varianee  of  zero 
will  eause  a  50  pereent  error  rate.  Henee,  we  should  not  expect  any  improvement  from 
the  maximum-likelihood  RAKE  receiver  designed  for  AWGN  and  pulse-noise  interfer¬ 
ence  when  hard  decision  Viterbi  decoding  is  used. 

4,  Soft  Decision  Viterbi  Decoding 

For  soft  decision  Viterbi  decoding,  can  be  written  as  [9] 


p‘{^-p)  PM 

i=o  V  i  J 


(2.76) 


where  p  is  the  duty  cycle  of  the  pulse-noise  interference,  and  P^(/)  is  the  probability  of 


selecting  a  code  word  a  Hamming  distance  d  from  the  correct  code  word,  given  that  i  of 
the  d  receptions  experience  interference.  To  calculate  P^{i)  we  must  integrate 
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(2.77) 


0(5|^)  = 


X 


d-i 


\k=\  y  V  k=\  J 

over  the  Tikonov  pdf  to  eliminate  the  dependenee  on  (j)  as  shown  in  (2.58),  numerieally 
invert  the  Laplaee  transform,  and  integrate  the  result  from  0  to  oo .  As  ean  be  seen  from 
(2.77),  this  must  be  done  d  +  l  times.  We  will  use  the  same  rate  14,  eonstraint  length  9 
eonvolutional  eode  where  the  first  three  terms  of  (2.71)  are  337^2  +2817^4  +21797^g , 
whieh  are  suffieient  to  ealculate  the  probability  of  bit  error.  Even  though  the  terms  are 


inereasing,  the  7^  terms  in  (2.71)  go  to  zero  mueh  faster.  An  exponential  MIP  is  assumed 

where  eaeh  linger  has  10  dB  less  signal  energy  than  the  previous  finger.  Figure  2.10 
shows  results  obtained  for  a  three  finger  RAKE  with  soft  deeision  Viterbi  deeoding  where 
^bl^Q- 15  dB,  the  pilot  tone  has  a  20-dB  loop  SNR,  a  ratio  of  direet-to-diffuse  power  of 
10  dB,  and  a  interferenee  signal  that  has  a  one  pereent  duty  eyele  {p  =  0.01 ),  a  ten  per- 
eent  duty  eyele  ( /?  =  0. 1 ),  and  a  100  pereent  duty  eyele  (/?  =  !).  The  ease  where  the  inter- 
ferenee  signal  is  on  all  the  time  was  also  simulated  for  low  values  of  Ef^jN j  sinee  the 

union  bound  given  by  (2.71)  is  know  to  be  loose  for  low  SNR.  Reeall  that  the  signal  input 
to  the  RAKE  reeeiver  is  assumed  to  be  a  Gaussian  random  variable  with  a  Rieean  distrib¬ 
uted  mean.  Multiplieation  by  a  Rieean  distributed  tap  weight  estimate  produees  a  noneen- 
tral  ehi-squared  mean  value  [21].  The  Matlab  statisties  toolbox  has  a  noneentral  ehi- 
squared  random  variable  generator.  These  samples  must  be  multiplied  by  cos(2^z))  where 
is  a  Tikonov  distributed  random  variable.  A  Tikonov  random  variable  generator  was 
developed  speeially  for  this  simulation.  The  mean  value  is  then  multiplied  by  -1  to  simu¬ 
late  all-zeros  being  transmitted.  This  eliminates  the  need  to  eonvolutionally  eneode  the 
data  sinee  an  all  zero  file  input  to  a  eonvolutional  eneoder  will  produee  an  all-zero  output 
file.  The  varianee  of  the  reeeived  Gaussian  random  variable  is  cr^  +  cr^  .  Sinee  the  mean 

and  varianee  are  known,  the  reeeived  Gaussian  random  variable  ean  now  be  simulated 
using  the  Matlab  randn  eommand.  Eaeh  finger  was  simulated  separately  and  summed 
together.  The  resultant  signal  was  then  quantized  and  eonvolutionally  deeoded  using  the 
Matlab  eommunieations  toolbox.  Any  data  bit  ones  at  the  deeoder  output  represent  a  bit 
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error,  so  the  bit  error  rate  was  eomputed  by  summing  the  values  in  the  eonvolutional  de- 
eoder  output  file  and  then  dividing  by  the  total  number  of  bits  simulated.  Monte-Carlo 
simulations  were  eondueted  for  eaeh  point  where  an  asterisk  is  shown  in  Figure  2.10.  A 
line  between  simulated  values  has  been  added  to  help  the  eye  follow  the  progression  be¬ 
tween  simulated  points.  Sinee  the  error  rate  varied  slightly  between  simulations,  three 
simulation  were  run  at  eaeh  point,  and  the  median  value  was  plotted  in  Figure  2.10.  Due 
to  the  eomputational  eomplexity,  only  error  rates  with  analytieal  values  above  10^^  were 
simulated.  These  simulation  results  are  found  to  be  within  1  dB  of  the  analytieal  results, 
whieh  shows  that  using  only  the  first  three  terms  of  the  union  bound  provides  an  aeeurate 
estimate  of  the  system  performanee. 

We  see  that  for  soft  deeision  Viterbi  deeoding,  pulsing  the  interferenee  signal  de¬ 
grades  the  performanee  mueh  more  than  it  does  for  hard  deeision  Viterbi  deeoding.  Soft 
deeision  Viterbi  deeoding  uses  soft  deeision  reeeiver  outputs  with  eonfidenee  estimates 
and,  if  any  ineorreet  deeisions  have  a  large  eonfidenee  value,  then  soft  deeision  deeoder 
performanee  signifieantly  degrades.  Essentially,  the  soft  decision  Viterbi  decoding  algo¬ 
rithm  attempts  to  combine  several  soft  decision  outputs  together  in  order  to  decode  the 
received  data.  The  maximum-likelihood  receiver  weights  the  soft  decision  outputs  by  the 
inverse  of  the  variance,  which  we  derived  in  Section  F,  to  avoid  placing  high  confidence 
values  on  bits  with  interference.  Unfortunately,  conventional  RAKE  receivers  do  not  per¬ 
form  this  function.  By  modeling  the  soft  decision  receiver  output  as  a  Gaussian  random 
variable  and  noting  that  the  sum  of  independent  Gaussian  random  variables  is  a  Gaussian 
random  variable  where  the  mean  is  the  sum  of  the  means  and  the  variance  is  the  sum  of 
the  variances  [15],  it  is  easy  to  see  that  any  bits  that  experience  interference  will  dramati¬ 
cally  increase  the  bit  error  ratio  unless  the  signal  power  is  greater  than  the  interference 
power. 

For  the  case  of  the  maximum-likelihood  RAKE  receiver  designed  for  AWGN  and 
pulse-noise  interference,  the  mean  and  variance  of  the  signal  on  each  RAKE  finger  both 
approach  zero  when  a  large  interference  signal  is  present.  This  means  that  soft  decision 
Viterbi  decoding  will  perform  much  better  in  the  presence  of  pulse-noise  interference 
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Figure  2.10  Coherent  RAKE  Receiver  with  Soft  Decision  Viterbi  Decoding. 

since  this  RAKE  receiver  does  not  allow  the  bits  with  interference  to  dominate  the  decod¬ 
ing  decision.  Results  for  a  three-finger  maximum-likelihood  RAKE  receiver  designed  for 
AWGN  and  pulse-noise  interference  with  soft  decision  Viterbi  decoding  where 
EblN^  =  15  dB,  the  pilot  tone  has  a  20-dB  loop  SNR,  a  ratio  of  direct-to-diffuse  power  of 

10  dB,  and  a  interference  signal  with  a  one  percent  duty  cycle  (/?  =  O.Ol) ,  a  ten  percent 

duty  cycle  (/:>  =  0.  l)  and  a  1 00  percent  duty  cycle  (/?  =  1 .0)  are  shown  in  Figure  2.11. 

As  expected.  Figure  2.1 1  shows  that  the  maximum-likelihood  RAKE  receiver  designed 
for  AWGN  and  pulse-noise  interference  mitigates  the  effect  of  pulse-noise  interference. 
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Figure  2.11  Maximum-Likelihood  RAKE  Receiver  for  AWGN  and  Pulse-Noise 
Interference  with  Soft  decision  Viterbi  decoding. 


When  the  interference  signal  is  on  ah  the  time,  the  performance  of  the  maximum- 
likelihood  RAKE  receiver  designed  for  AWGN  and  pulse-noise  interference  is  the  same 
as  the  conventional  RAKE  receiver. 

5,  Optimal  Finger  Spacing 

Since  we  assumed  an  exponential  multipath  intensity  profile  where  the  SNR  for 
each  finger  was  10  dB  below  the  previous  finger,  it  may  be  possible  to  obtain  better  per¬ 
formance  by  trying  to  reduce  the  difference  in  the  signal  power  on  each  finger.  Diversity 
performance  improves  when  we  decrease  the  finger  spacing  to  reduce  the  difference  in 
finger  SNR;  however,  if  we  limit  the  RAKE  to  three  fingers  then  the  total  received  power 
decreases  as  we  decrease  the  finger  spacing.  In  Figure  2.12  we  show  an  example  of  the 
multipath  intensity  profile  as  a  function  of  delay  with  two  different  possible  RAKE  finger 

spacing.  In  the  one  case  the  first  finger  collects  90  percent  of  the  received  power,  the  sec- 
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ond  nine  percent,  and  the  third  0.9  percent;  while  in  the  other  case,  the  first  finger  collects 
53.54  percent  of  the  received  power,  the  second  collects  24.87  percent,  and  the  third  fin¬ 
ger  collects  1 1.56  percent.  The  second  case  provides  better  diversity  performance  since 
the  received  power  on  the  three  fingers  are  closer  together,  but  it  only  collects  about  90 
percent  of  the  received  signal  power  as  opposed  to  the  first  case  which  collects  99.9  per¬ 
cent  of  the  received  signal  power. 


Figure  2.12  Multipath  Intensity  Profile  with  two  Different  RAKE  Finger  Spacings. 


In  order  to  find  the  optimal  tradeoff  between  losing  received  signal  power  and  getting 
better  performance,  we  used  hard  decision  Viterbi  decoding  and  started  with  a  RAKE  re¬ 
ceiver  that  recovered  99.9  percent  of  the  received  signal  power  and  compared  the  results 
with  a  RAKE  receiver  that  recovered  one  percent  less  power.  We  continued  reducing  the 
received  signal  power  by  one  percent  until  the  probability  of  bit  error  started  to  increase. 
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We  found  that,  for  a  three  finger  RAKE,  recovering  92  percent  of  the  received  signal 
power  produces  the  lowest  probability  of  bit  error.  In  Figure  2.13  we  show  the  results  for 
a  RAKE  receiver  that  recovers  99.9  percent  of  the  received  power  and  a  RAKE  receiver 
that  recovers  92  percent  of  the  received  power  for  hard  decision  Viterbi  decoding  where 
EblN^  =  15  dB,  the  pilot  tone  has  a  20-dB  loop  SNR,  a  ratio  of  direct-to-diffuse  power  of 


Figure  2.13  Two  Different  Finger  Spacings  with  Flard  Decision  Viterbi  Decoding. 


10  dB,  a  interference  signal  that  has  a  one  percent  duty  cycle  (/?  =  0.0 1) ,  ten  percent 
duty  cycle  (/?  =  O.l) ,  and  a  100  percent  duty  cycle  (/?  =  l.O) .  We  see  that  for  this  E^INq 
there  is  very  little  improvement  at  low  values  of  E^jN  ■  ,  and  the  improvement  at  high 
values  of  E^jN j  will  probably  not  be  noticed  in  a  practical  system.  The  hard  decision 
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results,  however,  show  that  reeovering  92  pereent  of  the  reeeived  power  produees  the 
lowest  probability  of  bit  error  for  large  E^jN j  . 

Sinee  the  maximum-likelihood  RAKE  reeeiver  for  AWGN  and  pulse-noise  inter- 
ferenee  with  soft  deeision  Viterbi  deeoding  reduees  the  probability  of  bit  error  for  pulse- 
noise  interferenee  down  to  the  irredueible  error  floor  eaused  by  phase  noise,  we  eompare 
the  performanee  of  a  finger  spaeing  that  reeovers  99.9  pereent  of  the  reeeived  power 
against  a  finger  spaeing  that  reeovers  92  pereent  of  the  reeeived  power  with  a  interferenee 
signal  that  is  on  all  the  time.  We  make  the  same  assumptions  that  we  did  for  hard  deei- 
sion  Viterbi  deeoding;  however,  with  soft  deeision  Viterbi  deeoding  the  error  rate  slightly 
inereases  when  only  92%  of  the  reeeived  power  is  reeovered  as  shown  in  Figure  2.14. 

It  seems  that  when  pulse-noise  interferenee  is  present,  performanee  with  both  hard 
and  soft  deeision  Viterbi  deeoding  does  not  improve  mueh  when  the  RAKE  finger  spae- 
ing  is  deereased  to  provide  better  diversity  performanee.  In  faet,  the  soft  deeision  per- 
formanee  was  slightly  worse.  These  results  were  obtained  with  =  15  dB.  The  re¬ 
sults  might  be  different  for  larger  or  smaller  ;  however,  based  on  the  results  ob¬ 

tained,  we  ean  eonelude  that  it  is  not  worth  the  extra  system  eomplexity  to  eonstantly 
measure  the  multipath  intensity  profile  and  adjust  the  finger  spaeing  for  large  Ef^jNj  . 


6.  Increased  Phase  Noise 

The  previous  results  have  all  been  for  a  loop  SNR  of  15  dB.  For  a  first  order 
phase-loeked  loop,  the  loop  SNR  ean  be  written 


Loop  SNR 


(2.78) 


where  E^  is  the  energy  of  the  pilot  tone,  cr^  is  the  loop  noise  power,  whieh  ineludes  both 
thermal  noise  and  interferenee  noise,  and  is  the  loop  bandwidth.  Deereasing  the  en- 
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Figure  2.14  Two  Different  Finger  Spacings  with  Soft  Decision  Viterbi  Decoding. 


ergy  of  the  pilot  tone  or  increasing  either  the  loop  bandwidth  or  noise  power  will  de¬ 
crease  Loopsjjp^ .  Figure  2.15  shows  the  significant  error  rate  increase  for  a  maximum- 

likelihood  RAKE  receiver  designed  for  AWGN  and  pulse-noise  interference  with  soft  de¬ 
cision  Viterbi  decoding  with  a  loop  SNR  of  10  dB,  15  dB  and  20  dB  when 
^bl^o  =  15  dB,  the  ratio  of  direct-to-diffuse  power  is  10  dB,  and  the  interference  signal 
has  a  ten  percent  duty  cycle.  Monte-Carlo  simulations  performed  using  MATLAB  show 
the  error  rate  is  0.5  for  a  loop  SNR  of  10  dB,  and  is  almost  10^^  at  15-dB  loop  SNR. 

Phase  noise  causes  an  irreducible  error  floor  and  problems  with  the  union  bound  since 
Pi2  ~  Pi4  ~  Pi6  ■  The  simulation  results  with  10-dB  and  15-dB  loop  SNR  have  lower  error 
rates  than  the  results  obtained  analytically  since  the  union  bound  used  to  obtain  analytic 
results  is  very  loose  when  phase  noise  dominates  system  performance.  Still,  the  simulated 
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Figure  2.15  Soft  Decision  Viterbi  Decoding  with  Various  Loop  SNR  Values. 


bit  error  ratio  (BER)  was  extremely  poor  for  a  decoded  signal.  Since  soft  decision  Viterbi 
decoding  cannot  provide  reliable  communication  for  a  loop  SNR  of  15  dB  or  below,  the 
classic  interference  strategy  of  interference  the  pilot  tone  is  highly  effective  in  this  case. 

Unlike  soft  decision  Viterbi  decoding,  hard  decision  Viterbi  decoding  only  uses 
the  receiver  estimate  of  what  bit  was  transmitted  without  any  estimate  of  the  receiver 
confidence.  While  phase  noise  increases  the  bit  error  rate,  a  convolutional  code  can  per¬ 
form  extremely  well  when  the  encoded  error  rate  is  less  than  10^^ .  Soft  decision  Viterbi 
decoding,  however,  uses  an  estimate  of  the  receiver  confidence,  which  is  adversely  af¬ 
fected  by  phase  noise.  Flence,  it  is  expected  that  hard  decision  Viterbi  decoding  will  per¬ 
form  better  for  10-dB  and  15-dB  loop  SNR.  Figure  2.16  shows  the  performance  for  a 
maximum-likelihood  RAKE  receiver  designed  for  AWGN  and  pulse-noise  interference 
with  hard  decision  Viterbi  decoding,  the  direct-to-diffuse  power  is  10  dB,  and  the  inter- 
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ference  signal  has  a  ten  percent  duty  cycle.  The  results  show  that  reliable  communication 
is  possible  even  with  10-dB  loop  SNR  when  hard  decision  Viterbi  decoding  is  used.  In 
fact,  the  performance  for  10-dB  and  15-dB  loop  SNR  was  better  for  hard  decision  Viterbi 
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Figure  2.16  Hard  Decision  Viterbi  Decoding  with  Various  Loop  SNR  Values. 


decoding  for  all  than  for  soft  decision  Viterbi  decoding.  When  the  loop  SNR  is 

equal  to  20  dB,  however,  soft  decision  Viterbi  decoding  performs  better  than  hard  deci¬ 
sion  Viterbi  decoding  for  small  EhjNj  .  The  soft  decision  Viterbi  decoding  used  to  obtain 

the  results  plotted  in  Figure  2.15  assumed  a  maximum-likelihood  receiver  for 
AWGN  and  pulse-noise  interference.  If  a  conventional  RAKE  receiver  designed  for 
AWGN  only  is  used,  then  the  performance  for  soft  Viterbi  decision  decoding  for  small 
EijjN-  is  significantly  worse  than  for  hard  decision  for  a  20-dB  loop  SNR.  A  decrease  in 
the  duty  cycle  of  the  pulse-noise  interference  signal  causes  problems  for  the  RAKE  re- 
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ceiver  optimized  for  AWGN  when  soft  deeision  Viterbi  deeoding  is  used;  however,  when 
hard  deeision  Viterbi  deeoding  is  used,  the  performanee  for  small  Ef^jNj  improves  as  the 

duty  eyele  deereases. 

The  maximum-likelihood  RAKE  reeeiver  for  pulse-noise  interferenee  was  de¬ 
rived,  and  this  deteetor  improves  the  performanee  for  soft  deeision  Viterbi  deeoding 
down  to  the  phase  noise  error  floor.  Unfortunately,  soft  deeision  Viterbi  deeoding  per- 
formanee  eannot  be  improved  below  the  phase  noise  error  floor.  Hard  deeision  Viterbi 
deeoding,  on  the  other  hand,  performs  mueh  better  than  soft  deeision  Viterbi  deeoding  for 
low  values  of  loop  SNR  sinee  soft  deeision  Viterbi  deeoding  suffers  more  from  the  ef- 
feets  of  imperfect  phase  synchronization.  A  plot  of  the  phase  error  as  a  function  of  loop 
SNR  was  shown  in  Figure  2.3.  As  the  loop  SNR  decreases,  the  phase  error  increases, 
which  leads  to  performance  degradation  as  a  result  of  imperfect  phase  synchronization. 
Since  soft  decision  Viterbi  decoding  makes  bit  decisions  by  combining  a  sequence  of  re¬ 
ceived  bits,  random  losses  due  to  imperfect  phase  synchronization  affect  the  entire  se¬ 
quence  of  bits  and  are  devastating  when  the  pilot  tone  has  a  loop  SNR  of  15  dB  or  less. 
Hard  decision  Viterbi  decoding  makes  decisions  on  a  bit-by-bit  basis  and  is  only  affected 
by  phase  noise  if  the  imperfect  phase  synchronization  loss  is  large  enough  to  cause  a  bit 
error. 

7.  Varying  the  Value  of  K 

The  previous  results  were  for  a  direct-to-diffuse  power  ratio  of  ten  (A  =  10) . 

There  are  other  important  cases  such  as  when  K  ^co  which  represents  no  fading  and 
K  =  0  which  represents  Rayleigh  fading.  The  RAKE  receiver  is  not  optimal  for  channels 
without  fading  because  all  the  signal  power  will  be  on  the  first  finger  and  the  other  fin¬ 
gers  will  only  add  additional  noise.  Worst  case  fading  occurs  when  there  is  no  direct  line 
of  sight  component  ( A  =  O) . 

In  Figure  2.17,  we  show  the  performance  for  a  three  finger  RAKE  receiver  with 
soft  decision  Viterbi  decoding  where  Ef^jN^  =  15  dB,  the  pilot  tone  has  a  20-dB  loop 

SNR,  the  interference  signal  has  a  100  percent  duty  cycle  (/?  =  l.O) ,  and  ratio  of  direct- 
to-diffuse  power  of  0,  10  and  100.  The  performance  obtained  for  a  direct-to-diffuse 
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power  ratio  of  100  approaches  that  obtained  for  K  ,  and  performance  is  only 
slightly  better  than  the  if  =  10  moderate  fading  case.  Worst  case  fading  if  =  0  signifi¬ 
cantly  increases  the  probability  of  bit  error.  Similar  results  are  found  when  the  interfer¬ 
ence  signal  is  pulsed  with  a  10  percent  duty  cycle  /?  =  0.1  as  shown  in  Figure  2.18. 


Figure  2.17  iL^ifE  Receiver  Performance  for  p  =  \.0  and  K=0, 10,  and  100. 

Figures  2.17  and  2.18  show  that  there  is  a  diminishing  return  gained  when  the 
value  of  K  increases  by  a  factor  of  10  from  10  to  100.  As  the  ratio  of  K  increases,  the 
power  on  the  first  finger  increases  and  decreases  on  the  second  and  third  fingers.  In  order 
to  compare  the  7f  =  100  with  results  from  an  AWGN  only  channel,  we  must  convert  the 

results  as  a  function  of  EtjNj  values  into  one  of  +  A^o)  •  Since  the  calculations 

were  performed  for  A’^/A^o  =15dB,  Ej^jNj  =4  dB  converts  to  +  Ag  j  =  3.65  dB. 
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Figure  2.18  RAKE  Receiver  Performance  for  /?  =  0.1  and  K=0,10  and  100. 


The  probability  of  bit  error  for  £'*/A^o  =  3-65  dB  AWGN  only  channel  and  a  constraint 
length  9  convolutional  code  in  [22]  is  slightly  lower  than  the  value  shown  at 
Ei^jNj  =  4  dB  in  Figure  2.17.  Since  the  RAKE  receiver  is  not  optimal  for  a  nonfading 

AWGN  channel  and  performance  for  A  =  100  approaches  performance  for  A  ^  oo  ,  the 
results  presented  in  this  dissertation  show  good  agreement  with  previously  published  re¬ 
sults. 
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III.  NONCOHERENT  RAKE  RECEIVER  PERFORMANCE 


The  performance  of  a  coherent  RAKE  receiver  degrades  when  the  pilot  tone  is  re¬ 
ceived  with  a  low  SNR,  so  we  now  turn  our  attention  to  the  noncoherent  RAKE  receiver 
which  does  not  require  any  knowledge  of  the  received  signal  phase.  In  this  dissertation, 
only  noncoherent  BFSK  is  considered.  The  transmitted  signal  for  a  data  bit  “0”  must  be 
orthogonal  to  the  transmitted  signal  for  a  data  bit  “1”  just  as  for  a  conventional  noncoher¬ 
ent  signal.  However,  the  transmitted  signal  must  also  be  orthogonal  to  delayed  versions 
of  the  transmitted  signal  when  the  delay  is  greater  than  the  RAKE  linger  spacing  for  a 
noncoherent  RAKE  receiver. 

For  the  noncoherent  RAKE  receiver,  the  transmitted  signal  for  a  bit  “0”  is  defined 
as 

a(t)  =  4lA^c  +  (j)) ,  (3.1) 

and  the  transmitted  signal  for  a  bit  “1”  is 

a(l)  =  42A,c(t)  cos((Ujt  +  (5'),  (3.2) 

where  42 A^  is  the  unfaded  carrier  amplitude,  c(t)  is  the  chipping  waveform,  which  has 
an  autocorrelation  function  defined  in  (2.5),  and  (j)  and  d  are  random  starting  phases. 

The  transmitted  frequencies  co^  and  a>^  are  chosen  to  be  orthogonal  over  the  bit  duration 
T  so  that: 


»T 

cos((»ot)cos((»jt)  =  0. 


(3.3) 


When  this  signal  is  transmitted  over  the  frequency-selective  fading  channel  shown  in 
Figure  1 .4,  the  received  signal  can  be  expressed  as 

(0  =  Z  V24«,c  (t  -  iT, )  cos  [a,  {t  -  iT^ )].  (3.4) 

i=i 

Each  of  the  /  components  is  modeled  as  a  flat  fading  signal  and  can  have  a  different  fad¬ 
ing  coefficients  a. .  When  a  bit  “0”  is  transmitted  a>^  =  when  a  bit  “1”  is  trans¬ 


mitted  .  The  received  signal  can  be  recovered  by  the  noncoherent  RAKE  receiver 

shown  in  Figure  3.1 
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Figure  3 . 1  Noncoherent  RAKE  Receiver. 

The  locally  generated  signals  are  5*0  (t)  =  V2c(t)cos((yot)  and  S';  (t)  =  V2c(t)cos((yjt) . 
At  the  output  of  the  mixer  on  the  linger  we  have 

i=\ 

for  the  received  signal  mixed  with  and 
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y2{^)  =  Y4'^A(^f{t-iT^)cos[(o^{t-iT^)y{t-jT^)cos[Q},{t-jT^)\,  (3.6) 

(=1 

for  the  received  signal  mixed  with  iSj  (t) .  The  integrator  then  integrates  these  signals, 
which  will  eliminate  any  delayed  components  of  the  received  signal  that  do  not  have  the 
same  delay  as  the  locally  generated  chipping  waveform  c(t),  since  the  chipping  wave¬ 
form  has  ideally  zero  correlation  with  delayed  versions  of  itself  as  shown  in  (2.5).  In  real¬ 
ity,  there  are  no  PN  sequences  that  give  zero  correlation  with  delayed  versions  of  them¬ 
selves,  but  practical  PN  sequences  can  be  long  enough  that  the  correlation  is  zero  for  all 
practical  purposes.  The  integrator  also  removes  the  double  frequency  term.  Since  the 
transmitted  signal  is  orthogonal,  in  the  absence  of  noise  the  output  of  the  integrator  on  the 
signal  branch  will  be  A^a.  and  the  non-signal  branch  integrator  output  will  be  zero.  The 

integrator  outputs  are  squared  and  combined.  The  results  of  the  correlation  of  the  re¬ 
ceived  signal  with  (t)  are  then  subtracted  from  the  corresponding  correlation  with 

iSg  (t) ,  and  the  receiver  decides  that  a  bit  “0”  was  transmitted  if  the  result  is  less  than 
zero.  Otherwise,  the  receiver  decides  a  bit  “1”  was  transmitted. 

A,  NONCOHERENT  RAKE  RECEIVER  ANALYSIS  WITH  PULSE-NOISE 
INTERFERENCE  AND  RICEAN  FADING 

We  will  now  derive  the  performance  of  a  noncoherent  RAKE  receiver  for  a  signal 
transmitted  over  a  Ricean  fading  channel  with  pulse-noise  interference.  Both  the  RAKE 
receiver  optimized  for  AWGN  channels,  analyzed  in  this  section,  and  the  maximum- 
likelihood  RAKE  receiver  optimized  for  pulse-noise  interference  ,  analyzed  in  the  next 
section,  will  be  investigated.  Convolutional  coding  with  both  hard  and  soft  decision 
Viterbi  decoding  will  also  be  considered. 

We  assume  without  loss  of  generality  that  the  signal  corresponding  to  bit  “0”  was 
transmitted.  The  conditional  probability  density  function  of  the  random  variable  at 

the  output  of  the  quadratic  detector  for  the  finger  of  the  signal  branch  given  a  signal 
amplitude  is  [17] 
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^\k  T 


where  1^  (•)  represents  the  modified  Bessel  function  of  the  first  kind  and  order  zero.  The 

average  received  power  on  the  k‘^  finger  is  al .  The  subscript  \k  is  used  to  identify  the 

k‘^  finger  of  the  signal  branch,  while  the  finger  of  the  noise  only  branch  will  be  la¬ 
beled  with  the  subscript  2k .  The  fading  channel  is  modeled  by  assuming  to  be  a  Ri- 

cean  random  variable.  The  probability  density  function  of  the  Ricean  random  variable  is 
[17] 


where  is  the  average  power  of  the  direct  signal  component  and  2Q^  is  the  average 
power  of  the  diffuse  signal  component  on  the  k‘^  finger.  The  total  received  power  on  the 
k‘''  finger  is  al  +  2Q^ .  The  conditioning  of  on  is  removed  by  integrating  over 


the  pdf  for  : 


fx,,  [Xxk )  =  I /x„  {^Xk  I  a, )  (a, )  da,. 


Substituting  (3.7)  and  (3.8)  into  (3.9),  we  obtain 


fx,,  (-^u)  “ 


:  ja,  exp 


Qz  ,  2  2 

k+O^k^k 

2alQl 


2alnl+alcrl 


2cTlkll 


J  J 

Jq  2 

.  ^  ^  ; 


(3.10) 


The  above  expression  can  be  evaluated  using  the  relationship 

4(x)  =  -^o(7-x), 

where  j  =  ^/^  ,  and  (6.633.4)  from  [18],  which  gives 


(3.11) 
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00 

J exp ( -aa,  )lo{/3a,)Jo{ra,)da,  =  —Qxp 


T  f^l 
\2a]' 


(3.12) 


Evaluating  (3.10)  using  (3.11)  and  (3.12),  we  get 


fx.^  {^\k) 


2a:,+al  2(2nJ+<T^)l  “(2n;+cr;J' 


(3.13) 


The  pdf  at  the  output  of  the  summer  is  then  given  by  the  eonvolution  of  the  pdf 
from  all  the  fingers  of  the  RAKE  receiver.  In  order  to  avoid  doing  these  convolutions  di¬ 
rectly,  the  Laplace  transforms  of  the  pdf  from  each  finger  are  multiplied  together  and  the 
inverse  transform  of  this  product  yields  the  pdf  at  the  output  of  the  summer.  The  Laplace 
transform  of  the  pdf  for  is  obtained  from 


(■^)  =  J /x„  {^ik ) exp {-sx,, ) dx,,. 


(3.14) 


Substituting  (3.13)  into  (3.14),  we  get 


2(2Ql+alf''\2(2Ql+al) 


By  making  the  change  of  variables 


T  = 


(3.15) 


(3.16) 


(3.17) 


we  can  evaluate  the  integral  in  (3.15)  using  (6.631.4)  from  [18]; 

[x"^^  QXY>{-ax^\j  {Px)dx  =  — ^  exp  -  — 

{  ^  ’  (2«f‘  t  4« 


(3.18) 
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Hence,  the  Laplace  transform  of  the  pdf  for  the  k‘^  finger  of  the  signal  branch  can  be 
written  as 


l  +  2s(2nl+al) 


exp 


-2s  a 


(i+2^(2n;+CTj)) 


(3.19) 


and  the  Laplace  transform  at  the  output  of  the  summer  for  the  signal  branch  is  simply  the 
product  of  the  k  signal  branch  fingers 

=  (3.20) 

k=\ 


The  pdf  for  the  noise  only  branch  can  be  found  from  (3.19)  by  setting  al  =kl\=0: 


2a: 


-exp 


X. 


2k 


J 


(3.21) 


As  previously  mentioned,  the  subscript  2k  identifies  the  k‘''  finger  on  the  noise  only 
branch,  i.e.,  a^  =  al  when  only  AWGN  is  present  and  a]  =al+  a^j  when  pulse-noise 
interference  is  present.  Using  (3.14),  we  obtain  the  Laplace  transform  of  (3.26)  as 


Fn  (*)  = 


1 


2a^s  + 1 


(3.22) 


The  pdf  at  the  output  of  the  summer  on  the  noise  only  branch  is  the  inverse  Laplace  trans¬ 
form  of  the  product  of  all  Laplace  transforms  for  all  the  fingers  at  the  input  of  the  sum¬ 
mer  as 

=  (3.23) 

k=l 

1.  Soft  Decision  Viterbi  Decoding 

The  probability  of  bit  error  for  a  system  with  convolutional  coding  is  upper 
bounded  by  [9] 

00 

P,<'Lb,P„  (3.24) 

d=dfree 
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where  is  the  total  number  of  nonzero  information  bits  in  all  weight  d  paths  and  is 
the  probability  of  seleeting  a  eode  sequenee  that  is  a  Hamming  distanee  d  away  from  the 
eorreet  eode  sequenee.  The  interleaver  randomizes  the  order  of  the  bits  with  interferenee 
at  the  deeoder  input,  so  ean  be  written  [9]  as 


/=0 


(d 


v'y 


p‘{i-prpM 


(3.25) 


where  p  is  the  duty  eyele  of  the  pulse  noise  interferenee,  and  P^  (i)  is  the  probability  of 

seleeting  a  eode  sequenee  a  Hamming  distanee  d  from  the  eorreet  eode  sequenee  given 
that  i  of  the  d  reeeptions  experienee  interferenee.  For  a  rate  14,  eonstraint  length  9  eonvo- 
lutional  eode,  the  first  three  nonzero  terms  of  the  series  are  33/^2  +281/^4  +2179/^g , 

whieh  are  suffieient  to  aceurately  eompute  /),  when  P,<  10  ^ .  To  ealeulate  Pj{i)  when 


d-i  bits  are  reeeived  with  only  AWGN  and  i  bits  experienee  interferenee,  we  must 
evaluate 


^  (0  “  (^1  )  \  f2x{,^2y^^2 

0  1_0 

In  order  to  evaluate  (3.26),  we  obtain  the  required  pdfs  from 

f,Xx,)  =  L-^[F,Xs)] 


dXy 


and 


where 


and 


(^) = I  )  1  [  tipx,,  I  = ^o+^j) 

V  k=i 


y-v  L 

y  V  k=i 


^ = f  n  ^  I  ^ = ^0 )  1  [  n  ^  I  = ^0 + ) 

V  lc=l 


r  L 

n 

y  V  k=\ 


I) 

'I 


(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 
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The  evaluation  of  (3.26)  requires  numerieal  inversion  of  (3.29)  and  (3.30),  followed  by 
numerieal  integration  of  (3.26).  The  inner  integral  does  not  have  to  be  evaluated  numeri- 
eally  sinee,  as  is  well  known  [23], 


(3.31) 


The  required  inversions  are  done  numerieally  using  the  method  discussed  in  [24],  which 
inverts  the  Laplace  transform  O  (5)  at  any  value  of  x  from 


/W 


exp(ax) 

2x 


Red)(a)  +  2^(-l)”0 

n=\ 


r 

nni 

Q  + 

V 

X 

)_ 

(3.32) 


where  a  is  a  number  greater  than  the  real  part  of  any  singularity  of  d)(^) .  The  trapezoidal 
rule  is  then  used  to  integrate  the  numerically  inverted  data  points  to  obtain  Pci{i)  ■  These 
numerical  inversions  must  be  done  d  +  l  times.  For  example,  the  /[j  includes  the 

case  where  no  soft  decision  outputs  experience  interference,  where  one  experiences  inter¬ 
ference,  where  two  experience  interference,  and  so  on  to  the  case  where  all  twelve  out¬ 
puts  experience  interference.  These  results  are  then  substituted  into  (3.25)  to  obtain  . 

Figure  3.2  shows  the  probability  of  bit  error  for  a  three  RAKE  optimized  for 
the  AWGN  channel  with  rate  14,  constraint  length  9,  Viterbi  soft  decision  Viterbi  decod¬ 
ing  when  =  15  dB,  the  ratio  of  direct  to  diffuse  power  is  10  dB  and  the  interference 

signal  has  a  one  percent  duty  cycle  (/?  =  O.Ol) ,  ten  percent  duty  cycle  (/?  =  0.  l) ,  and  a 
interference  signal  that  is  on  all  the  time  (/?  =  1 .0) .  Results  show  that  the  probability  of 
bit  error  increases  when  the  interference  duty  cycle  p  is  decreased.  This  is  due  to  the  fact 
that  the  interference  power  increases  as  the  duty  cycle  decreases,  since  =  Nj  IrT,^ 
where  N j  is  the  interference  noise  power  spectral  density  and  7],  is  the  bit  duration.  Soft 

decision  Viterbi  decoding  uses  a  sequence  of  soft  decision  receiver  outputs  to  make  a  bit 
decision.  If  any  of  these  soft  decision  receiver  outputs  have  a  large  variance,  then  the 
probability  of  bit  error  will  significantly  increase. 
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Figure  3.2  Noncoherent  RAKE  Receiver  Optimized  for  the  AWGN  Channel  Perform¬ 
ance  with  Rate  'A,  Constraint  Length  9,  Viterbi  Soft  Decision  Viterbi  Decoding  as  a 

Function  of  Duty  Cycle  p. 

The  bits  with  interference  will  also  affect  future  bit  decisions,  which  will  cause  additional 
decoding  errors.  Comparing  Figure  3.2  with  Figure  2.10,  shows  that  the  coherent  RAKE 
receiver  performs  better  for  low  values  of  E^jN .  but  that  the  noncoherent  RAKE  re¬ 
ceiver  performs  below  the  irreducible  error  floor  caused  by  phase  noise  in  the  coherent 
RAKE  receiver  at  high  values  of  E^^jN .  .  The  noncoherent  combining  losses  cause  the 

noncoherent  RAKE  receiver  to  have  catastrophic  performance  for  E^jNj  below  5  dB 
when  the  interference  signal  is  on  all  the  time  and  approximately  for  Ef^jNj  below  26  dB 
when  the  interference  signal  had  a  one  percent  duty  cycle  {p  =  0.01) .  It  should  also  be 
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noted  that  pulsing  the  interference  signal  increased  the  probability  of  bit  error  for  both  the 
coherent  and  noncoherent  RAKE  receivers  designed  for  AWGN  channels.  Therefore,  we 
can  conclude  that  pulse-noise  interference  is  a  very  effective  interference  strategy  against 
a  RAKE  receiver  designed  for  AWGN  channels  with  soft  decision  Viterbi  decoding.  As  a 
result,  pulse-noise  interference  will  be  very  effective  against  commercial  systems  since 
all  commercial  RAKE  receivers  are  designed  in  this  manner. 

2,  Hard  Decision  Decoding 

If  hard  decision  Viterbi  decoding  is  utilized,  in  (3.24)  is  [9] 


d  (d 

P,-  z  \  . 


;_£/+! 


v) 


(3.33) 


when  d  is  odd,  and 


P  -1  2 

d 


■P^  -(1-^)2 


d  (d 


(3.34) 


+  Z  .  -p  ■[}-p) 


when  d  is  even.  In  (3.33)  and  (3.34),  the  channel  transition  probability  p  is 


p[Pb\^^  +  P)[Ph\^^  =^l)^  (3-35) 

where  p  is  the  duty  cycle  of  the  pulse-noise  interference,  (^/),  |  cr^  =  erg +cr]^  is  the 
probability  of  bit  error  when  a  single  bit  experiences  interference,  and  (^/),  |  cr^  =  cTq  j  is 

the  probability  of  bit  error  when  a  single  bit  is  transmitted  with  only  AWGN  present. 

Figure  3.3  shows  the  probability  of  bit  error  for  a  three-finger  RAKE  with  rate  14 
constraint  length  9  Viterbi  hard  decision  decoding  when  E^jN^  =  15  dB,  the  ratio  of  di¬ 
rect  to  diffuse  power  is  10  dB  and  the  interference  signal  has  a  one  percent  duty  cycle 
(/?  =  O.Ol) ,  ten  percent  duty  cycle  (/?  =  0.  l) ,  and  a  interference  signal  that  is  on  all  the 

time  (/?  =  1 .0) .  Soft  decision  Viterbi  decoding  uses  a  sequence  of  soft  decision  receiver 
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Figure  3.3  Noncoherent  RAKE  Receiver  with  Rate  V2  Constraint  Length  9  Viterbi 

Hard  Decision  Decoding. 


outputs  to  make  bit  decisions.  Note  that  interference  one  of  these  soft  decision  receiver 
outputs  may  cause  the  current  bit  decision  to  be  incorrect  as  well  as  several  of  the  follow¬ 
ing  bits  using  a  previous  soft  decision  receiver  output  with  interference.  Hard  decision 
Viterbi  decoding  makes  hard  decisions  at  the  receiver  and  does  not  provide  an  estimate  of 
the  confidence  in  the  bit  decision.  As  a  result,  hard  decision  Viterbi  decoding  treats  all  bit 
errors  equally  and  does  not  allow  a  single  bit  to  with  interference  dominate  the  decoding 
decision.  Further  note  that  the  hard  decision  receiver  outputs  have  few  errors  when  the 
pulse-noise  interference  duty  cycle  decreases.  Hence,  hard  decision  Viterbi  decoding  can 
provide  bit  error  probabilities  of  10 or  lower  when  the  pulse-noise  interference  duty 
cycle  is  one  percent,  which  is  a  significant  improvement  over  the  soft  decision  Viterbi 
decoding  that  required  Ej^jNj  »30  dB  to  achieve  this  bit  error  probability. 
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B,  MAXIMUM-LIKELIHOOD  RECEIVER  FOR  NONCOHERENT 
CHANNELS  WITH  GAUSSIAN  NOISE,  PULSE-NOISE 
INTERFERENCE,  AND  SOFT  DECISION  VITERBI  DECODING 

The  maximum-likelihood  receiver  performance  for  a  noncoherent  RAKE  receiver 
with  AWGN,  pulse-noise  interference,  and  soft  decision  Viterbi  decoding  is  derived  in 
this  section.  Since  in  (3.24)  is  equivalent  to  a  d‘’'  order  diversity  for  soft  decision 
Viterbi  decoding  [25],  the  analysis  in  this  section  can  be  applied  to  a  wide  range  of  sys¬ 
tems  that  use  diversity  as  well  as  error  correction  coding.  The  performance  of  the  coher¬ 
ent  RAKE  receiver  analyzed  previously  dramatically  improved  when  it  was  designed  for 
maximum-likelihood  performance  with  AWGN  and  pulse-noise  interference.  As  previ¬ 
ously  mentioned,  the  pulse-noise  interference  considered  will  have  a  duty  cycle  p  and 
interference  power  spectral  density  K  .j p  when  the  interference  signal  is  on.  Since  the 
interference  signal  is  modeled  as  additive  Gaussian  noise,  the  noise  variance  is  equal  to 
cTq  when  there  is  only  AWGN  present  and  <j\  -l-  crj  when  the  interference  is  on.  Assum¬ 
ing  that  the  diversity  receptions  are  independent,  a  reasonable  assumption  for  a  properly 
designed  system,  the  joint  density  function  is  the  product  of  the  marginal  density  func¬ 
tions  given  in  (3.7).  We  obtain  the  joint  pdf  for  the  “0”  bit  when  i  diversity  receptions 
experience  interference  and  d-i  diversity  receptions  have  only  Gaussian  noise  as 

(3.36) 

where  (t)  is  a  sequence  of  m  diversity  receptions  of  a  bit  “0”  and  a]  is  the  average 
received  power  for  the  c’'  RAKE  linger.  The  joint  pdf  for  a  bit  “1”  is  given  by 

(3.37) 
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(2  2  \  / 
[(J^  +  (Jj  )  m=d-M 
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60 


where  (tj  is  a  sequenee  of  m  diversity  receptions  for  a  bit  “1”.  The  likelihood  ratio  is 
given  by 


(3.38) 


Substituting  (3.36)  and  (3.37)  into  (3.38),  we  get 


f  j  V-‘ 
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exp 


2a, 


0  m=l 


d-i 


n^. 


m-\ 


exp 


-1 


u 


_ 1 

X 

'  2,2 

J  m=d-i+l  \ 

^  CTo+cr^.  j 

f  1  A 


d-i 


So  V  2^0  y 


exp 


^  Z-fo„(0+2«: 


2cr, 


0  m=l 
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a. 


(0 


(3.39) 
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J  m=d-i+\  \ 
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Taking  the  natural  log  and  rearranging  terms,  we  get 


d-i 


(0 
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(3.40) 
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m=d-i+l 
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2  2 

V  / 


Note  that  on  the  Ti  («)  branch  will  be  the  same  as  on  the  To  (^)  branch  since  it 
was  assumed  that  each  RAKE  linger  will  have  flat  fading.  Hence,  the  ln[/Q  (•)]  term 
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adds  the  same  value  on  both  sides  of  (3.40),  and  we  ean  further  simplify  the  maximum- 
likelihood  reeeiver  to 


cr. 


1  d-i 


m=l 


d  Si  1  d-i 

5o 


^1\  VV<  2  VV  ■  /  2 

j  m=d-i+\  m=\ 


4 


2  2 

^ j  I  m=d-i+l 


whieh  shows  that  the  output  of  eaeh  quadratie  deteetor  should  be  weighted  by  the  inverse 
of  the  reeeived  bit  varianee  +  aj  or  cTq  depending  on  whether  the  bit  experienees  in- 
terferenee  or  not,  respeetively.  A  diagram  of  the  maximum-likelihood  eombiner  for 
AWGN  and  pulse-  noise  interferenee  is  shown  in  Figure  3.4,  where  l/ cr^  =  l/ cTq  when 

only  AWGN  is  present  and  ij <j^  =  i/^cTq  +  <j^j  ^  when  the  reeeived  bit  experienees  inter¬ 
ferenee.  A  maximum-likelihood  noncoherent  RAKE  receiver  designed  for  AWGN  and 


Figure  3.4  Maximum-Likelihood  Combiner  for  AWGN  and  Pulse-Noise  Interfer¬ 
ence. 

pulse-noise  interference  is  shown  in  Figure  3.5.  Equation  (3.41)  shows  that  the  maxi- 
mum-likelihood  noncoherent  RAKE  receiver  scales  the  received  bit  by  the  inverse  of  the 
received  bit  variance  <jI  +  cr^  or  <jI  depending  on  whether  the  bit  experiences  interfer¬ 
ence  or  not,  respectively.  We  first  examine  the  effect  of  this  scaling  on  the  signal  branch 
in  order  to  analyze  the  maximum-likelihood  noncoherent  RAKE  receiver.  We  define 
^\k  ~  ^\kl ’  which  is  the  random  variable  on  the  A:'*  finger  of  the  signal  branch. 
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Figure  3.5  Maximum-Likelihood  RAKE  Reeeiver  for  AWGN  and  Pulse-Noise  Inter 

ference. 


To  find  the  pdf  for  Zj^ ,  we  make  the  transformation  Zj^. 


=  /cr^ ,  and  from  (3.7)  we 


f  (z  u)-‘e;:nr 

2  I  2ctI  I  ”  CT, 


(3.42) 


where  Iq  (•)  represents  the  modified  Bessel  funetion  of  the  first  kind  and  order  zero, 
is  assumed  to  be  a  Rieean  random  variable  defined  in  (3.14),  <jI  =  <jI  when  only  AWGN 
is  present  and  <jI  =<jI+<j^j  when  pulse-noise  interferenee  is  present.  We  eliminate  the 
dependenee  on  by  substituting  (3.42)  and  (3.9)  into  (3.10),  leading  to 


\k)  2q2  P  2q2 


:  j  a,  exp 
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(3.43) 


The  integral  ean  be  evaluated  using  (3.11)  and  (6.633.4)  from  [23],  whieh  yields 


f  (z  )-  _2a^+z^ 

l[l^t+CTl)\\2Cl^+CTl 


(3.44) 
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Making  the  substitutions  =  —r  and  =  — — ,  we  get 


If  ^Xk+^Pk  I  j 

n 


(3.45) 


whieh  is  (5)  in  [26].  To  ealeulate  the  Laplaee  transform  of  (3.45)  we  substitute  (3.45)  into 
(3.14)  whieh  yields 


0  dz,,.{3A6) 

1  +  4 


By  making  the  transformation  of  variables  u  =  and  du  =  dz^i^jl^J^  ,  we  get 

Z7  1  Pk  ,  1  X.z]  r  , 
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which  can  be  evaluated  using  (3.11)  and  (6.631.4)  from  [23]  ; 

P 


I  Qxp(^-ax^  ^  ( y^x)  dx  = 


exp 


V  4«y 


(3.48) 


Performing  the  integration  and  simplifying,  we  get 
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2'^(l  +  4)  +  l 


exp 


a 


2.(l  +  4) 


(3.49) 


The  Laplace  transform  at  the  output  of  the  signal  braneh  summer  is  the  product  of  the  L 
signal  branch  fingers,  so  that 

^'.v,(4  =  n^'x„(4  (3-50) 

^=1 


The  pdf  for  the  noise-only  braneh  of  the  maximum-likelihood  RAKE  reeeiver  ean 
be  found  from  (3.44)  by  setting  the  reeeived  direet  signal  power  to  zero  [a^  =  Oj  and  the 

reeeived  diffuse  signal  power  to  zero  {iQ.^  =  0^ .  Performing  this  substitution,  we  get 


/z,,(^2j  =  Texp 


(3.51) 


where  represents  the  random  variable  on  the  k‘^  linger  of  the  noise-only  braneh.  The 

Laplace  transform  for  the  noise-only  branch  is  calculated  by  substituting  (3.51)  into 
(3.14),  which  yields 


(^) 


1 

2s  +  \' 


(3.52) 


Since  the  maximum-likelihood  RAKE  receiver  seales  the  received  signal  by  the  inverse  of 
the  varianee  <jI  +  cr^  or  <jI  depending  on  whether  the  bit  experienees  interferenee  or  not, 
respeetively,  the  Laplace  transform  of  the  pdf  for  the  noise-only  RAKE  lingers  does  not 
depend  on  the  reeeived  signal  varianee  (cr^ )  •  The  Laplace  transform  of  the  pdf  at  the 
output  of  the  summer  is  the  produet  of  the  L  noise-only  lingers: 
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(3.53) 


k=\ 

The  Laplace  transform  at  the  output  of  the  signal  branch  summer  and  the  noise- 
only  branch  summer  must  be  inverted  and  substituted  into  (3.26)  to  calculate  the  bit  error 
probability.  As  an  example,  we  use  a  three  finger  RAKE  with  EjN^^  =15  dB,  and  a  10- 

dB  ratio  of  direct-to-diffuse  power  to  calculate  the  performance  of  the  maximum- 
likelihood  noncoherent  RAKE  receiver  against  a  interference  signal  with  one  percent  duty 
cycle  (/?  =  O.Ol) ,  ten  percent  duty  cycle  (/?  =  0.  l) ,  and  a  interference  signal  that  is  on  all 

the  time  (/?  =  1 .0) .  Results  presented  in  Figure  3.6  show  that  the  maximum-likelihood 

noncoherent  RAKE  receiver  with  rate  14,  constraint  length  9  soft  decision  Viterbi  decod¬ 
ing  improves  performance  when  pulse-noise  interference  is  present.  In  fact,  the  worst 
case  interference  at  low  SNR  values  was  achieved  when  the  interference  was  on  all  the 
time,  which  indicates  that  pulse-noise  interference  is  not  effective  against  a  maximum- 
likelihood  receiver  designed  for  both  AWGN  and  pulse-noise  interference.  The  perform¬ 
ance  of  the  maximum-likelihood  noncoherent  RAKE  receiver  is  significantly  better  than 
the  performance  of  the  standard  noncoherent  RAKE  receiver,  which  has  catastrophic  per¬ 
formance  for  Ef^jN j  below  25  dB  when  the  interference  signal  has  a  one  percent  duty 

cycle.  The  maximum-likelihood  noncoherent  RAKE  receiver  has  a  probability  of  bit  error 
at  or  below  10^^°  for  all  values  of  Ej^jN  ■  for  this  case.  In  summary,  a  noncoherent  RAKE 

receiver  designed  for  AWGN  performs  worse  as  the  pulse-noise  interference  duty  cycle 
decreases,  whereas  the  opposite  is  true  for  the  maximum-likelihood  noncoherent  RAKE 
receiver. 

Comparing  the  maximum-likelihood  coherent  RAKE  receiver  performance  shown 
in  Figure  2.1 1  with  that  of  the  maximum-likelihood  noncoherent  RAKE  receiver,  we  see 
that  the  noncoherent  case  does  not  have  an  irreducible  phase  noise  error  floor  and  is  only 
limited  by  AWGN  noise.  When  the  interference  signal  is  on  all  the  time  or  has  a  ten  per¬ 
cent  duty  cycle,  the  maximum-likelihood  coherent  RAKE  receiver  performs  better  for  low 
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Figure  3.6  Maximum-Likelihood  Noncoherent  RAKE  Receiver  with  Rate  V2  Con¬ 
straint  Length  9  Viterbi  Soft  Decision  Viterbi  Decoding. 


values  of  EbjNj ;  however,  when  the  interference  signal  has  a  one  percent  duty  cycle  the 
maximum-likelihood  noncoherent  RAKE  receiver  performs  better  for  ah  values  of 
EblNj  because  the  maximum-likelihood  coherent  RAKE  receiver  has  an  irreducible 

phase  noise  error  floor.  In  both  cases  the  worst  performance  occurs  when  the  interference 
signal  is  on  ah  the  time.  In  conclusion,  the  maximum-likelihood  RAKE  designed  for 
AWGN  and  pulse-noise  interference  provides  dramatically  better  performance  in  the 
presence  of  pulse-noise  interference  than  a  commercial  RAKE  receiver  designed  for 
AWGN  channels. 
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IV.  CONCLUSIONS 


A,  RESEARCH  SUMMARY 

The  performance  of  a  direct  sequence  spread  spectrum  (DSSS)  signal  transmitted 
over  a  Ricean  frequency-selective,  slowly  fading  channel  with  pule-noise  interference, 
phase  noise  and  AWGN  has  been  analyzed.  Both  coherent  and  noncoherent  modulation 
with  an  exponential  multipath  intensity  profile  were  considered.  These  signals  were  com¬ 
pared  using  a  direct  to  diffuse  power  ratio  of  ten  (A  =  lO) ,  which  represents  moderate 

fading.  The  performance  of  conventional  RAKE  receiver  optimized  only  for  AWGN  was 
compared  with  the  maximum-likelihood  RAKE  receiver  optimized  for  pulse-noise  inter¬ 
ference  and  AWGN.  A  rate  14  ,  constraint  length  9  convolutional  code  was  also  consid¬ 
ered  with  both  hard  and  soft  decision  Viterbi  decoding. 

Analysis  of  the  coherent  RAKE  receiver  with  soft  decision  Viterbi  decoding  dem¬ 
onstrated  how  critically  important  it  is  to  have  a  strong  pilot  tone.  In  Figure  2.15  the  per¬ 
formance  for  10-dB  loop  SNR  and  15-dB  loop  SNR  does  not  provide  reliable  communi¬ 
cation.  In  fact,  the  performance  for  lO-dB  loop  SNR  produced  a  catastrophic  50  percent 
error  rate  and  15-dB  loop  SNR  produced  a  10^^  error  rate.  In  order  to  avoid  these  disas¬ 
trous  error  rates,  commercial  wireless  systems  dedicate  as  much  as  20  percent  of  the  total 
power  to  the  pilot  tone.  The  20-dB  loop  SNR  performance  had  an  irreducible  phase  noise 
error  floor  slightly  above  10  ,  which  is  more  than  sufficient  for  wireless  communica¬ 

tion  systems.  Fiber  optic  systems  which  operate  at  error  rates  of  approximately  10 
would  require  a  loop  SNR  greater  than  15  dB. 

Commercial  coherent  RAKE  receivers  designed  for  AWGN  with  soft  decision 
Viterbi  decoding  perform  poorly  in  the  presence  of  pulse-noise  interference.  Results  pre¬ 
sented  in  Figure  2.10  demonstrate  that  the  commercial  coherent  RAKE  receiver  performs 
worse  as  the  duty  cycle  of  the  pulse-noise  interference  signal  decreases.  A  soft  decision 
receiver  provides  estimates  of  whether  the  received  signal  is  a  bif’O”  or  a  bit  “1”  and  an 
estimate  of  the  confidence  in  this  decision.  When  a  strong  interference  signal  is  present. 
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the  soft  decision  receiver  can  produce  a  bit  error  with  a  high  confidence  of  being  the  cor¬ 
rect  estimate.  The  convolutional  decoder  will  produce  a  series  of  bit  errors  when  the 
minimum  weight  path  through  the  trellis  includes  a  bit  error  with  a  high  confidence  esti¬ 
mate.  Since  the  peak  interference  power  increases  as  the  pulse-noise  interference  duty 
cycle  decreases,  it  is  logical  to  expect  the  bit  error  rate  to  increase  as  the  duty  cycle  of  the 
pulse-noise  interference  duty  cycle  decreases.  The  maximum-likelihood  co\\qxqyARAKE 
receiver  optimized  for  pulse-noise  interference  and  AWGN  was  derived  and  the  analysis 
proved  that  measuring  the  variance  of  the  received  signal  on  a  bit-by-bit  basis  and  scaling 
each  bit  by  the  inverse  of  the  variance  is  optimal  according  to  the  maximum-likelihood 
criteria.  The  efficacy  of  this  receiver  was  demonstrated  by  the  fact  that  it  reduced  the 
probability  of  bit  error  down  to  the  irreducible  phase  noise  error  floor  when  pulse-noise 
interference  was  present.  Results  were  derived  analytically  up  to  the  point  where  no  ana¬ 
lytical  solution  to  the  required  inverse  Laplace  transform  existed,  so  numerical  techniques 
were  used  both  for  inversion  and  integration  to  calculate  the  probability  of  bit  error.  Ana¬ 
lytical  results  for  the  RAKE  receiver  designed  only  for  AWGN  were  verified  by  simula¬ 
tion  for  the  case  when  the  interference  signal  was  on  all  the  time  and  the  simulation  re¬ 
sults  were  within  1  dB  of  the  analytical  results.  The  maximum-likelihood  coherent  RAKE 
receiver  designed  for  AWGN  and  pulse-noise  interference  was  also  simulated  for  10-dB 
and  15-dB  loop  SNR.  The  simulated  probability  of  bit  error  for  10-dB  and  15-dB  loop 
SNR  were  below  the  analytic  values  because  the  union  bound  is  loose  when  phase  noise 
dominates  system  performance. 

Results  showed  that  the  RAKE  receiver  designed  for  AWGN  only  and  hard  deci¬ 
sion  Viterbi  decoding  performed  better  when  the  duty  cycle  of  the  interference  signal 
was  reduced.  A  hard  decision  receiver  only  provides  an  estimate  of  whether  the  received 
data  is  a  bit  “0”  or  a  bit  “1”  and  does  not  provide  an  estimate  of  the  received  bit  decision 
confidence.  As  a  result,  a  bit  that  experiences  interference  will  not  dominate  the  decoding 
decision  since  the  decoder  does  not  use  estimates  of  the  received  bit  confidence.  Decreas¬ 
ing  the  duty  cycle  of  the  interference  signal  also  decreases  the  number  of  bits  that  ex¬ 
perience  interference.  With  hard  decision  Viterbi  decoding,  fewer  bits  that  experience 
interference  at  the  input  to  the  decoder  produce  fewer  errors  at  the  output  of  the  decoder. 
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In  Figure  2.9,  when  the  duty  eyele  of  the  interferenee  signal  was  one  pereent,  the  prob¬ 
ability  of  bit  error  was  approximately  10 or  lower  for  all  values  of  E^^jN j  .  This  is 

signifieantly  lower  than  the  probability  of  bit  error  for  the  same  ease  with  soft  deeision 
Viterbi  deeoding,  whieh  produeed  a  eatastrophie  50  pereent  error  rate  when  Ej^jNj  was 

less  than  3  dB.  From  Figure  2.10,  soft  deeision  Viterbi  deeoding  required  Ef^jNj  =  15  dB 

to  reaeh  a  10 probability  of  bit  error.  Henee,  hard  deeision  Viterbi  deeoding  is  a  good 
method  to  limit  the  effeet  of  pulse-noise  interferenee. 

In  addition,  hard  deeision  Viterbi  deeoding  also  limits  the  effeets  of  phase  noise. 
The  results  showed  that  soft  deeision  Viterbi  deeoding  ean  not  provide  reliable  eommuni- 
eation  for  either  10-dB  or  15-dB  loop  SNR.  However,  hard  deeision  Viterbi  deeoding  ean 
reduee  the  probability  of  bit  error  to  10  for  25  dB  or  greater  E^jNj  with  a  loop  SNR  of 

10  dB,  and  the  15-dB  loop  SNR  ease  only  required  1 1  dB  Ef^jNj  .  This  shows  that  inter¬ 
fering  with  the  pilot  tone  would  not  eause  a  eatastrophie  error  rate  when  the  pilot  tone 
and  data  are  transmitted  on  different  frequeneies  and  hard  deeision  eonvolutional  eoding 
is  used. 


Changing  the  RAKE  linger  spaeing  in  order  to  improve  the  diversity  performanee 
by  redueing  the  differenee  in  the  signal  power  on  eaeh  linger  was  investigated.  Sinee  the 
analysis  assumed  a  three-finger  RAKE,  when  the  finger  spaeing  was  redueed  to  deerease 
the  differenee  in  finger  SNR,  the  total  reeeived  power  also  deereased.  Repeated  analysis 
for  hard  deeision  Viterbi  deeoding  demonstrated  that  reeovering  92  pereent  of  the  re¬ 
eeived  power  produeed  the  lowest  probability  of  bit  error.  Deereasing  the  finger  spaeing 
further  to  provide  better  diversity  performanee  did  not  fully  eompensate  for  the  additional 
reeeived  signal  power  reduetion.  When  soft  deeision  Viterbi  deeoding  was  analyzed  with 
a  RAKE  linger  spaeing  that  reeovered  92  pereent  of  the  reeeived  signal  power,  the  prob¬ 
ability  of  bit  error  was  approximately  the  same  as  the  RAKE  reeeiver  that  reeovered  99.9 
pereent  of  the  reeeived  signal  power  and  eaeh  finger  is  10  dB  below  the  previous  one.  It 
should  be  noted  that  this  analysis  is  speeifie  to  a  three-finger  RAKE  reeeiver  and  it  is  pos¬ 
sible  that  a  RAKE  reeeiver  with  four  fingers  or  more  eould  benefit  by  adapting  the  finger 
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spacing.  Additionally,  a  15-dB  Ej^jN^  was  assumed  in  the  analysis.  Adapting  the  RAKE 
finger  spacing  could  be  more  effective  for  higher  levels  of  AWGN. 

Noncoherent  RAKE  receivers  designed  for  AWGN  with  soft  decision  Viterbi  de¬ 
coding  perform  miserably  in  the  presenee  of  pulse-noise  interferenee.  Results  presented 
in  Figure  3.4  demonstrate  that  the  noneoherent  RAKE  receiver  performance  degrades  as 
the  duty  eycle  of  the  pulse-noise  interferenee  signal  deereases.  When  the  pulse-noise  in¬ 
terference  signal  uses  a  one  pereent  duty  eycle,  the  noneoherent  RAKE  reeeiver  with  soft 
deeision  Viterbi  deeoding  has  eatastrophie  performanee  for  values  of  E^jNj  below  26 

dB.  By  comparison,  at  26-dB  E^jN j  the  eoherent  RAKE  reeeiver  designed  for  AWGN 

had  a  probability  of  bit  error  slightly  above  10 for  a  interferenee  duty  cyele  of  one  per¬ 
eent  when  the  pilot  tone  was  reeeived  with  a  20-dB  loop  SNR.  A  strong  pilot  tone  is  eriti- 
cal  for  coherent  systems.  The  noncoherent  RAKE  reeeiver  performed  better  than  the  co¬ 
herent  RAKE  reeeiver  when  the  pilot  tone  was  received  with  a  10-dB  loop  SNR,  whieh 
had  a  50  percent  error  rate  for  all  values  of  E^jNj  .  Unlike  the  coherent  RAKE  receivers, 

the  noncoherent  RAKE  receiver  does  not  have  an  irreducible  probability  of  bit  error  floor. 
Since  noncoherent  receivers  demodulate  the  received  bits  without  using  a  pilot  tone,  they 
are  immune  to  random  phase  fluctuations  that  cause  problems  for  coherent  receivers. 

Hard  decision  Viterbi  deeoding  also  limits  the  effect  of  pulse-noise  interferenee 
on  a  noncoherent  RAKE  reeeiver  designed  for  AWGN.  When  soft  deeision  Viterbi  deeod¬ 
ing  is  used,  the  probability  of  bit  error  inereases  as  the  duty  cyele  of  the  interferenee  sig¬ 
nal  deereases;  however,  the  opposite  is  true  for  hard  deeision  Viterbi  deeoding.  Sinee  a 
hard  deeision  reeeiver  does  not  provide  an  estimate  of  the  eonfidence  of  the  bit  deeision, 
the  decoder  will  not  have  the  problem  of  an  erroneous  bit  decision  with  a  high  confidence 
estimate  of  being  correet.  The  eoherent  RAKE  reeeiver  designed  for  AWGN  with  hard 
deeision  Viterbi  decoding  and  20-dB  loop  SNR  performed  better  than  the  noncoherent 
RAKE  receiver  for  all  values  of  E^^jNj  .  A  strong  received  pilot  tone  is  neeessary  for  the 

coherent  RAKE  receiver  to  outperform  the  noneoherent  RAKE. 
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The  maximum-likelihood  noneoherent  RAKE  receiver  for  pulse-noise  interference 
and  AWGN  was  also  derived.  Like  the  coherent  case,  measuring  the  variance  of  the  re¬ 
ceived  signal  on  a  bit-by-bit  basis  and  scaling  each  bit  by  the  inverse  of  the  variance  is 
optimal  according  to  the  maximum-likelihood  criteria.  Our  analysis  demonstrated  that, 
when  the  interference  signal  has  a  10  percent  duty  cycle,  this  receiver  has  a  probability  of 
bit  error  below  10 for  all  values  of  ,  and  for  a  one  percent  duty  cycle  the  prob¬ 
ability  of  bit  error  was  below  10^^'’  for  all  values  of  .  Since  the  maximum- 

likelihood  noncoherent  RAKE  receiver  for  pulse-noise  interference  and  AWGN  effec¬ 
tively  mitigates  pulse-noise  interference,  the  interference  signal  would  have  to  transmit 
continuously;  however,  for  values  of  E^jNj  above  10  dB  the  probability  of  bit  error  is 

below  10^^ .  Therefore,  the  maximum-likelihood  noncoherent  RAKE  receiver  for  pulse- 
noise  interference  and  AWGN  is  highly  resistant  to  interference. 

In  summary,  a  RAKE  receiver  designed  for  AWGN  does  not  perform  well  in  the 
presence  of  pulse-noise  interference.  The  maximum-likelihood  RAKE  receiver  designed 
for  AWGN  and  pulse  noise  interference  mitigates  the  effect  of  pulse-noise  interference 
for  both  coherent  and  noncoherent  RAKE  receivers.  This  demonstrates  how  important  it 
is  to  design  the  receiver  for  the  intended  operational  environment.  It  is  also  critically  im¬ 
portant  for  a  coherent  RAKE  receiver  to  have  a  strong  pilot  tone.  Random  fluctuations  of 
the  pilot  tone  caused  by  high  levels  of  AWGN  have  a  devastating  effect  on  soft  decision 
Viterbi  decoding.  Hard  decision  Viterbi  decoding  can  limit  the  impact  of  phase  noise,  but 
soft  decision  Viterbi  decoding  performs  better  than  hard  decision  Viterbi  decoding  when 
the  pilot  tone  is  received  with  high  SNR. 

B,  SUGGESTIONS  FOR  FURTHER  RESEARCH 

It  was  recently  shown  that  the  multipath  amplitude  distribution  for  indoor  chan¬ 
nels  is  well  modeled  by  the  Nakagami-m  distribution  [27].  Since  direct  sequence  spread 
spectrum  with  binary  phase-shift  keying  or  quadrature  phase-shift  keying  is  often  used 
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for  indoor  channels,  it  would  be  interesting  to  analyze  the  eoherent  and  noncoherent 
RAKE  reeeiver  performanee  when  the  signal  is  transmitted  over  a  Nakagami  fading 
ehannel. 

This  dissertation  assumed  a  rate  V2  eonvolutional  eode  that  is  one  of  two  convolu¬ 
tional  codes  used  in  IS95  eellular  eommunication  systems.  For  a  noisy  reverse  ehannel, 
IS95  systems  use  a  rate  1/3  code,  while  the  rate  14  convolutional  code  considered  in  this 
dissertation  is  specified  for  the  forward  ehannel.  The  interested  reader  ean  repeat  this  re¬ 
search  for  a  rate  1/3,  constraint  length  9  convolutional  code.  The  new  3G  standards  spec¬ 
ify  both  a  convolutional  code  and  a  turbo  eode.  It  would  be  very  intriguing  to  compare 
eoherent  RAKE  reeeivers  with  eonvolutional  eoding  and  turbo  coding  for  an  AWGN 
channel  with  pulse-noise  interferenee. 

The  noncoherent  RAKE  reeeiver  in  this  dissertation  used  binary  frequency-shift 
keying  (FSK).  In  [14]  it  was  shown  that  a  differential  phase  shift  keyed  (DPSK)  signal 
can  be  represented  as  an  orthogonal  signal  over  a  period  of  two  bits.  Thus  DPSK  typi- 
eally  has  a  3-dB  advantage  over  FSK  for  noncoherent  eommunications.  The  noneoherent 
RAKE  receiver,  however,  must  be  uneorrelated  from  finger  to  finger,  and  in  [9]  it  is 
shown  that  a  noneoherent  RAKE  receiver  using  DPSK  modulation  would  require  two 
ehips  between  fingers  instead  of  the  one  ehip  between  fingers  that  FSK  requires.  It  is  not 
intuitively  obvious  that  using  DPSK  modulation  instead  of  FSK  modulation  would  im¬ 
prove  the  performanee  of  the  noneoherent  RAKE  reeeiver,  sinee  the  DPSK  modulation 
would  require  the  noneoherent  RAKE  reeeiver  to  operate  at  twice  the  ehipping  rate  of  the 
noncoherent  RAKE  reeeiver  with  FSK  modulation  for  the  same  ehannel. 
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APPENDIX  A 


The  Gauss-Chebyshev  quadrature  method  given  in  (2.59)  ean  be  used  as  a  quiek 
iterative  teehnique  to  ealeulate  the  Q  funetion.  Mathematieally,  the  Q  funetion  ean  be  ex- 

” r  1  ^ 

pressed  as  Q{x)  =  [  , —  e  ^  du.  This  integral  of  the  Gaussian  pdf  is  used  to  determine 

the  probability  of  bit  error  for  BPSK  signals  transmitted  over  an  AWGN  ehannel.  While 
no  elosed  form  solution  exists,  the  Q  funetion  has  been  tabulated  for  various  values  of  x. 
We  ean  ealeulate  the  value  for  Q(2)  by  substituting  -2  for  m  and  1  for  <j^  in  (2.46), 
whieh  is  the  Laplaee  transform  of  the  Gaussian  pdf  The  Gauss-Chebyshev  quadrature 
method  ean  now  be  used  to  invert  the  Laplaee  transform  and  integrate  it  from  zero  to  in¬ 
finity.  Sinee  we  ehose  -2  for  the  mean  value,  the  result  should  be  the  same  as  Q(2).  The 
following  MATLAB  eode  implements  the  Gauss-Chebyshev  quadrature  method; 


v=200; 

i=sqrt(-l); 

sum=0; 

0=2.5; 

for  k=l:(v/2) 

tk=tan((((2*k)-l)*pi)/(2*v)); 

sum=sum-l-real(LAP(o+(i*o*tk)))-l-(tk*imag(LAP(o+i*o*tk))); 

end 

res=sum/v; 


funotion[y]=LAP(s) 
y=exp(-2  *  s+(0 .5  *  s^2)); 

When  the  eode  is  run  the  result  is  0.0228,  whieh  is  the  same  value  given  for  Q(2)  in  [14]. 

In  this  dissertation,  the  Laplaee  transforms  that  were  inverted  required  ohoosing 
different  values  for  c  as  the  signal-to-noise  ratio  ohanged.  In  general,  small  values  of  c 
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were  used  for  low  signal-to-noise  ratios  and  larger  values  were  ehosen  for  higher  signal- 
to-noise  ratios.  The  following  MATLAB  eode  was  used  for  the  eoherent  RAKE  reeeiver 
optimized  for  AWGN  with  p  =  \: 


d=2*pi*besseli(0,3 1 .623); 

v=200; 

i=sqrt(-l); 

for  p=l:41 

ebnj=(10^(.l*(p-l))); 
nj=1.585/ebnj; 
sum=0; 
for  k=l:v/2 

tk=tan((2*k-l)*pi/(2*v)); 
integ=0; 
for  1=1:99 
o=.05; 
ifp>2 
0=0.1; 
end 
ifp>8 
0=0.5; 
end 
ifp>25 
0=1.88; 
end 

a(l)=EX(o+(i*o*tk),((l*0.01*pi)-pi),nj); 

end 

for  m=l:98 

g(m)=(a(m)+a(m+ 1  ))/2; 


76 


n=cos((((m+m+l)/2)*0.01*pi)-pi); 
z(m)=0.02*pi*g(m)*(exp(3 1 .623*n))/d; 
integ=integ+z(m); 
end 

sum=sum+real(integ)+tk*imag(integ); 

end 

res(p)=sum/v; 


funetion[y  ] =F(s  ,p ,  ej  ) 

No=  1; 

Nj=.l+((2*ej)/l); 

K=10; 

g=l/22; 

B=2*1.585*(eos(p)^2)*(l-((s/2)*No))*s; 

B2=2*1.585*(eos(p)^2)*(l-((s/2)*Nj))*s; 

y=((((exp(-2*K*0.9*B2*g/(l+2*0.9*B2*g)))/(l+2*g*0.9*B2))*((exp(- 

2*K*0.09*B2*g/(l+2*0.09*B2*g)))/(l+2*g*0.09*B2))*((exp(- 

2*K*0.009*B2*g/(l+2*0.009*B2*g)))/(l+2*0.009*B2*g)))^12)*((((exp(- 

2*K*0.9*B*g/(l+2*0.9*B*g)))/(l+2*g*0.9*B))*((exp(- 

2*K*0.09*B*g/(l+2*0.09*B*g)))/(l+2*g*0.09*B))*((exp(- 

2*K*0.009*B*g/(l+2*0.009*B*g)))/(l+2*0.009*B*g)))^0); 


This  ealeulates  the  term.  To  ealeulate  we  change  the  function  F  as  shown  in  the 
following: 


function[y]=F(s,p,ej) 
No=  1; 

Nj=l+((2'^ej)/l); 

K=10; 
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g=l/22; 

B=2*1.585*(cos(p)^2)*(l-((s/2)*No))*s; 

B2=2*1.585*(cos(p)^2)*(l-((s/2)*Nj))*s; 

y=((((exp(-2*K*0.9*B2*g/(l+2*0.9*B2*g)))/(l+2*g*0.9*B2))*((exp(- 

2*K*0.09*B2*g/(l+2*0.09*B2*g)))/(l+2*g*0.09*B2))*((exp(- 

2*K*0.009*B2*g/(l+2*0.009*B2*g)))/(l+2*0.009*B2*g)))^14)*((((exp( 

2*K*0.9*B*g/(l+2*0.9*B*g)))/(l+2*g*0.9*B))*((exp(- 

2*K*0.09*B*g/(l+2*0.09*B*g)))/(l+2*g*0.09*B))*((exp(- 

2*K*0.009*B*g/(l+2*0.009*B*g)))/(l+2*0.009*B*g)))^0); 


Next,  can  be  calculated  using  the  following  function  F: 


function[y]=F(s,p,ej) 

No=  1; 

Nj=.l+((2*ej)/l); 

K=10; 

g=l/22; 

B=2*1.585*(cos(p)^2)*(l-((s/2)*No))*s; 

B2=2*1.585*(cos(pr2)*(l-((s/2)*Nj))*s; 

y=((((exp(-2*K*0.9*B2*g/(l+2*0.9*B2*g)))/(l+2*g*0.9*B2))*((exp(- 

2*K*0.09*B2*g/(l+2*0.09*B2*g)))/(l+2*g*0.09*B2))*((exp(- 

2*K*0.009*B2*g/(l+2*0.009*B2*g)))/(l+2*0.009*B2*g)))^16)*((((exp( 

2*K*0.9*B*g/(l+2*0.9*B*g)))/(l+2*g*0.9*B))*((exp(- 

2*K*0.09*B*g/(l+2*0.09*B*g)))/(l+2*g*0.09*B))*((exp(- 

2*K*0.009*B*g/(l+2*0.009*B*g)))/(l+2*0.009*B*g)))^0); 


Finally,  the  probability  of  bit  error  can  be  calculated  from 
F,=(33/;,)  +  (281/;,)  +  (2179/;,). 
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APPENDIX  B 


The  numerical  Laplace  transform  inverse  given  in  (3.32)  can  be  implemented  as 
shown  below.  This  example  inverts  the  Laplace  transform  sj +1^  ,  which  has  a  time 

domain  inverse  (t/2)sin(t) .  When  t  =  njl ,  the  time  domain  function  has  a  value  of 

02111,  which  is  the  same  result  we  get  using  the  following  MATLAB  code  to  invert  the 
Laplace  transform  at  the  point  t  =  ;r/2  ; 


i=sqrt(-l); 

v=1000; 

t=pi/4; 

a=7/t; 

inside=0; 

for  k=l:v 

c=a+(k*pi*i/t); 
ifrem(k,2)  ==  0 

inside=inside+(real(w(c))); 

else 

inside=inside+(real((- 1 )  *  (w(c)))); 
end 
end 

tmp=real(w(a))+(2*  inside); 

tmp2=(exp(a*t)/(2*t)); 

tmp3=tmp2*tmp; 

function[y]=w(s) 

y=s/((s^2+l)^2); 


The  numerical  Laplace  transform  inversions  performed  for  this  dissertation,  required  the 
constant  a  to  be  1  Ht  instead  oil  It  as  shown  in  the  previous  example.  The  following  code 
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was  used  to  invert  the  Laplaee  transforms  for  the  noneoherent  RAKE  reeeiver  designed 
for  AWGN  with  rate  V2,  eonstraint  length  9  eonvolutional  coding  and  Viterbi  soft  deci¬ 
sion  decoding: 

i=sqrt(-l); 
v=2000; 
for  p=l:40 

ebnj=(10^(.l*(p-l))); 
nj=l+(31.623/(0.1*ebnj)); 
for  j= 1:200 
t=j*6.0*(nj/3 1.623); 
a=ll/t; 
inside=0; 
for  k=l:v 

c=a-l-(k*pi*i/t); 
if  rem(k,2)==0 

inside  =  inside+(real(w(c,nj))); 
else 

inside  =  inside+(real((-l)*w(c,nj))); 
end 
end 

tmp=real(  w(a,nj  ))-l-(2  *  inside) ; 
tmp2=(exp(a*t)/(2*t)); 
res(j  )=tmp2  *  tmp ; 
end; 

for  zz=  1:200 
if  res(zz)<0; 

res(zz)=0; 

end 

end 

sum=0; 
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for  m=l;199 

int  1  =(res(m)+res(m+ 1  ))/2; 
int2=intl  *6.0*(nj/3 1 .623); 
sum=sum+int2; 
end 

endres(p)=sum; 

end 

For  the  p  =  \  ease,  we  ean  now  use  the  following  for  the  function  w  to  calculate  the 
term: 

function[y]=w(s,ej) 
sigl=0.9*3 1.623*  1/22; 
sig2=0.09*3 1.623*  1/22; 
sig3=0.009*3 1.623*  1/22; 
alphal=0.9*3 1.623*20/22; 
alpha2=0.09*3 1 .623*20/22; 
alpha3=0.009*3 1 .623*20/22; 

A=(l./((2.*ej.*s)+l)).^12; 

B=(l./((4.*s)+l)).^0; 

C=((l./(2.*s.*(2.*sigl+ej))).*(exp((- 

2.*s.*alphal)./(l+(2.*s.*(2.*sigl+ej))))).*(l./(2.*s.*(2.*sig2+ej))).*(exp((- 

2.*s.*alpha2)./(l+(2.*s.*(2.*sig2+ej))))).*(l./(2.*s.*(2.*sig3+ej))).*(exp((- 

2.*s.*alpha3)./(l+(2.*s.*(2.*sig3+ej)))))).^0; 

D=((l./(l+2.*s.*(2.*sigl+l))).*(exp((- 

2.*s.*alphal)./(l+(2.*s.*(2.*sigl+l))))).*(l./(2.*s.*(2.*sig2+l))).*(exp((- 

2.*s.*alpha2)./(l+(2.*s.*(2.*sig2+l))))).*(l./(2.*s.*(2.*sig3+l))).*(exp((- 

2.*s.*alpha3)./(l+(2.*s.*(2.*sig3+l)))))).^12; 

y=(C.*D).*(l-((l./s).*(A.*B))); 

The  value  can  now  be  calculated  by  using  the  following  for  the  function  w: 
function[y]=w(s,ej) 
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sigl=0.9*3 1.623*  1/22; 
sig2=0.09*3 1.623*  1/22; 
sig3=0.009*3 1.623*  1/22; 
alphal=0.9*3 1.623*20/22; 
alpha2=0.09*3 1 .623*20/22; 
alpha3=0.009*3 1 .623*20/22; 

A=(l./((2.*ej.*s)+l)).^14; 

B=(l./((4.*s)+l)).^0; 

C=((l./(2.*s.*(2.*sigl+ej))).*(exp((- 

2.*s.*alphal)./(l+(2.*s.*(2.*sigl+ej))))).*(l./(2.*s.*(2.*sig2+ej))).*(exp((- 

2.*s.*alpha2)./(l+(2.*s.*(2.*sig2+ej))))).*(l./(2.*s.*(2.*sig3+ej))).*(exp((- 

2.*s.*alpha3)./(l+(2.*s.*(2.*sig3+ej)))))).^0; 

D=((l./(l+2.*s.*(2.*sigl+l))).*(exp((- 

2.*s.*alphal)./(l+(2.*s.*(2.*sigl+l))))).*(l./(2.*s.*(2.*sig2+l))).*(exp((- 

2.*s.*alpha2)./(l+(2.*s.*(2.*sig2+l))))).*(l./(2.*s.*(2.*sig3+l))).*(exp((- 

2.*s.*alpha3)./(l+(2.*s.*(2.*sig3+l)))))).^14; 

y=(C.*D).*(l-((l./s).*(A.*B))); 

Finally,  can  be  obtained  using: 

funetion[y]=w(s,ej) 
sigl=0.9*3 1.623*  1/22; 
sig2=0.09*3 1.623*  1/22; 
sig3=0.009*3 1.623*  1/22; 
alphal=0.9*3 1.623*20/22; 
alpha2=0.09*3 1 .623*20/22; 
alpha3=0.009*3 1 .623*20/22; 

A=(l./((2.*ej.*s)+l)).^16; 

B=(l./((4.*s)+l)).^0; 

C=((l./(2.*s.*(2.*sigl+ej))).*(exp((- 

2.*s.*alphal)./(l+(2.*s.*(2.*sigl+ej))))).*(l./(2.*s.*(2.*sig2+ej))).*(exp((- 
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2.*s.*alpha2)./(l+(2.*s.*(2.*sig2+ej))))).*(l./(2.*s.*(2.*sig3+ej))).*(exp((- 

2.*s.*alpha3)./(l+(2.*s.*(2.*sig3+ej)))))).^0; 

D=((l./(l+2.*s.*(2.*sigl+l))).*(exp((- 

2.*s.*alphal)./(l+(2.*s.*(2.*sigl+l))))).*(l./(2.*s.*(2.*sig2+l))).*(exp((- 

2.*s.*alpha2)./(l+(2.*s.*(2.*sig2+l))))).*(l./(2.*s.*(2.*sig3+l))).*(exp((- 

2.*s.*alpha3)./(l+(2.*s.*(2.*sig3+l)))))).^16; 

y=(C.*D).*(l-((l./s).*(A.*B))); 

Lastly,  the  probability  of  bit  error  can  be  calculated  from 
/;=(33/;i)  +  (281/;,)  +  (2179/;,). 
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